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PREFACE 


The seven sections in this report summarize recent prediction 
techniques and results obtained by U.S. propagation reaarchers and system 
engineers- The results are written for use. by satellite! system engineers 
working fn the 10 to 100 GHz frequency range. 

.The key subjects covered in these sections are: 

• Attenuation, Section II 

• Path Diversity, Section III 

• Signal Fluctuations and Low Angle Fading,, Section IV 

• Depolarization Effects, Section V 

• Bandwidth Coherence, Section VI ■ 

• Sky Noise, Section VII 

Each section may be utilized individually and example calculations have been 
included as a guide for the system engineer. 
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1.0 PURPOSE 


i ■ 

i 

This report provides a compilation of prediction, methods and related 
propagation studies for the evaluation of earth-space paths operating above 
10 GHz. The topics covered in this report are: 

Attenuation 

Path Diversity 

SignasT Fluctuations and Low Angle Fading 

Depolarization Effects 

Bandwidth Coherence . 

Sky Noise 

The techniques described here have been developed from recent 
ongoing NASA supported studies and from the relevant published literature. 
These techniques represent the state of knowledge -ad verse .effects of 

the earth's- atmosphere, on reliable earth-space transmissions above 10 GHZ. 

■ - - o! ve: s Df‘- 

...This report provides propagation data in a format suitable for use 

by earth-space link system designers operating in the frequency range from 10 
to 100 GHz in the United States and Canada. In this frequency range the 
troposphere can have a significant .effect on the carrier to noise ratio of a 
propagating wave.. Typically, the troposphere attenuates- and depolarizes the 
carrier signal and adds broadband amplitude and phase noise to the signal. 

The resulting carrier-to-noise ratio reduction reduces the allowable data 
rate for a given bit error rate (digital systems) and the quality to 
transmission (analog systems). In the most severe cases- the medium will 
■significantly attenuate the carrier and destroy the transmission capabilities 
of the link (termed a link outage). The frequency of occurrence and average 
outage time per year are usually of most interest;to -system designers. 
Propagation studies to date now allow the predictions to be made with a high 
degree of certainty and have developed means to reduce the frequency and 
length of these outages. 
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ZQ ORGANIZATION OF THIS REPORT 


< This report is arranged in seven relatively independent sections 
covering the key topics related to the interaction of the troposphere and 
earth-space propagation paths. Each section attempts to present a 
description af the effect(s) and an example calculation related to a typical 
- communication system. Because these examples provide .a concise guide to the 
calculations required, an index to the examples is gi-yen in Table 1. 

- • An attempt has been made to keep the style are! format 

'-consistent,, however, this has not always been possible Three authors were 
involved and this document has evolved over a 6 month-jieriod. Thus, some 
'inconsistencies are inevitable and surely some arithmetic mistakes are still 
present. ’ 
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TABLE 1 


GUIDE TO PROPAGATION EXAMPLES 


Subject * ? tkrcation of Example 

Clear Air Attenuation Sact. II, Par. 2.0 

Rain Attenuation .Statistics (Global, Rice-Holmberg, Lin, 

Dutton-Qougherty model s ) SSct. II, Para. 3.5.2. 1 


Rain Attenuation .Statistics (given various 
statistfcal inputs) . . . ... . 

Path Diversity „ ....... 

Amplitude Fluctuations 

Phase and Angle-af-Arrival Variations . . 

Average Signal Gain Reduction ...... 

•Depolarization Statistics ........ 

Amplitude. Coherence., (rain induced) . . . 

Phase Coherence (ionospheric induced) . , 

- * * ' * ‘ 
Sky ; No1se Statistics. . . • 


Sect':. ill, Para. 3.5.3 
, . and 3.5.4 

-S&ct,. Ill, Para. 2.2 

Sect;' IV, Para. 6.1 

Sect .. IV, Para. 6.2 

3&ct: ‘IV, Para. 6.3 

.. V, Para. 3.0 

3&ct. VI, Para. 2.1.1 

. SSfct. VI, Para. 3.1 

, S&ct. VII, Para. 2.1.2 
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3..0 FREQUENCY BANOS FOR EARTH-SPACE 'COMPLICATION 

i 

i' 


Within the guidelines established by the Iritstmat-ional 
Telecommunications Union (ITU) for Region 2, the Federal Communications 
Commission. (FCC) in the U.S. and the Department of Conmuni cat ions (DOC) in 
Canada regulate the earth-space' frequency allocations- In most cases, the 
FCC and; DOC regulations are more restrictive than the ETU regulations for 
Region. 2. Therefore, the FCC and DOC earth-space- '.fregaency bands are 
considered in this report. 

The services which operate via an earth— spaoe link are listed in 
Table- .'2-. along; with- their ITU index number. The definitions of each of these 
services are giv.en in Ref. 1. The specific frequency allocations for these 
services are. listed in Table 3. These -frequency allocations are taken from 
the FCC’s Eight Notice of Inquiry and represent an estimate of' the Region 2 
frequency allocations beyond 1981 (.two years following the 1979 World 
Administrative Radio Conference). The Canadian freqisricy proposals are 
similar and are not included here. Both .frequency proposals are given in 
detail in References 2 and 3. 
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TABLE 2 


TELECOMMUNICATION SERVICES UTILIZING, 

i 

EARTH- SPACE PROPAGATION LINKS ■ 

Fixed Satellite (84AG) 

' Mobile Satellite (military satellite)" (8AAG8) 
Aeronautical Mobile Satellite (8AAE&); 
Maritime Mobile Satellite (84AGD) 

Land Mobile Satellite (84AGD) 

Broadcasting Satellite (84AP) 

1 kadionavigation Satellite (84AQ) 

: Earth Exploration Satellite (84ASAj) 
Meteorological Satellite (84.AT) 

Amateur Satellite (84ATA) 

Standard Frequency Satellite (.84AT8) 

Space Research (84ATD) 

Space Operation (84ATE) 

t 

Radio Astronomy (75A) 
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TABLE 3 


PROPOSED UNITED STATES, EARTH -SPACF FREQUENCY ALLOCATIONS 


FREQUENCY (GHz ) 


10.6 

- 10.7 

10.7 

- 12,2 

11.7 

- 12.5 

12.0 

- 13.25 

12.75 

- 14.2 

13.4 

- 14.0 

14.0 

- 14.5 

14.4 

- 14.5 

14.5 

- 15.35 

15.2 

- 15.14 

15.35 

- 15.4 

16.6 

- 17.1 

— « j * 

17.7 

- 17.9 

17.7 

- 21.2 

20.2 

- 21.2 

21.2 

- 21.4 

22.21 

- 22.5 


SERVICES 

RADIO ASTRONOMY 

FIXED SATELLITE (Downlink) 

BROADCASTING SATELLITE 
FIXED SATELLITE (Uplink) 

SPACE RESEARCH (Downlink! 

EARTH-EXPLORATION SATELLITE (Active Sensors) 
•STANDARD FREQUENCY SATELLITE (Uplink). 

FIXED SATELLITE (Uplink) 

SPACE .RESEARCH (Space-Space) 

SPACE RESEARCH 

EARTH-EXPLORATION SATELLITE (Passive) 

•SPACE RESEARCH (Passive) 

earth Exploration satellite (Active sensors) 
iPh fit HjfiiANCH (Aeiivfe Sensors) 
METEOROLOGICAL SATELLITE (Downlink) 

FIXED SATELLITE (Downlink)' ' ' ' 

MOBILE SATELLITE (Downlink) 

STANDARD FREQUENCY SATELLITE (Downlink) 

EARTH EXPLORATION SATELLITE (Passive) 

EARTH EXPLORATION SATELLITE (Passive) 

SPACE RESEARCH (Passive) 



TABLE 3 

FREQUENCY (GHz) 

22.6 23.6 

23.6 * 24 


24.0 - 24.05 
25.25 - 27.5 

25.25 - 26.25 
27.5 - 31.0 

30.0 - 31.0 

30 - 31.3 

31 - 31.3 

31.3 - 31.8 


3$: = am 

34.2 - 35.2 

36.0 - 37.0 

40.0 - 41.0 


(Continued) 

SE RVICE 

RADIO ASTRONOMY 
INTER-SATELLITE 

EARTH EXPLORATION SATELLITE (Passive) 

SPACE RESEARCH (Passive) 

RADIO ASTRONOMY 
AMATEUR SATELLITE 

STANDARD’ FREQUENCY SATELLITE .(Downlink) 

EARTH EXPLORATION SATELLITE ..(Space-to-Space) 
INTER- SATELLITE 
, FIXED SATELLITE (Uplink) 

MOBILE SATELLITE (Uplink) 

STANDARD FREQUENCY SATELLITE- (Downlink) 

SPACE RESEARCH 

SPACE RESEARCH (Passive) 

• EARTH EXPLORATION SATELLITE (Passive) 

' RADIO ASTRONOMY 

spAGE : rSesIeaHgH 

SPACE RESEARCH (Active) 

EARTH EXPLORATION SATELLITE (Active)’ 

SPACE RESEARCH (Passive) 

EARTH EXPLORATION SATELLITE (Passive) 

FIXED SATELLITE (Downlink) 

MOBILE SATELLITE (Downlink) 



TABLE 3 (Continued) 



49.8 

- 50.0 

50.0 

-51.4 

50,2 

- 50.4 

50.4 

- 51.4 

51.0 

- 51.4 

51.4 

- 59.0 

54^25 

- 53.2 

59.0 

64.0 

64.0 

- 66.0 

54.0 

- 65.0 

65.0 

- 66.0 

66.0 

- 71.0 


SERVICE; ' 

BROADCASTING SATELLITE 
FIXED SATELLITE (Uplink) 

MOBILE SATELLITE (Uplink) 

AERONAUTICAL MOBILE SATELLITE 
MARITIME MOBILE SATELLITE 
AERONAUTICAL RADIONAVIGATION SATELLITE 
MARITIME RADIONAVIGATION SATELLITE 
AMATEUR SATELLITE 
FIXED SATELLITE (Uplink) 

SPACE RESEARCH (Passive) 

EARTH - EXPLORATION SATELLITE (Passive) • 
MOBILE SATELLITE (Upllns) 

FIXED SATELLITE (Uplink) 

"EARTH EXPLORATION SATELLITE (Passive) 
SPACE RESEARCH (Passive) * - — 

iNIffWPLMfF 

INTERSATELLITE 

EARTH EXPLORATION SATELLITE (Passive) 
SPACE RESEARCH (Passive) 

SPACE RESEARCH 

AERONAUTICAL MOBILE SATELLITE 
MARITIME MOBILE SATELLITE 



FREQUENCY (GHz) 


TABLE 3 (Continued) 

SERVICE 


AERONAUTICAL RADIONAVIGATION SATELLITE 
MARITIME RADIONAVIGATION SATELLITE 


71.0 - 74,0 

FIXED SATELLITE (Uplink)' 
MOBILE SATELLITE (Uplink) 

76.0 - 77.0 

SPACE RESEARCH (Active) 

EARTH EXPLORATION SATELLITE (Active) 

81.0 - 84.0 

MOBILE SATELLITE (Downlink)' 

> 

FIXED SATELLITE (Downlink) 

84.0 - 86.0 

BROADCASTING SATELLITE 

86.0 - 92.0 

EARTH EXPLORATION SATELLITE (Passive) 
RADIO ASTRONOMY 


SPACE RESEARCH (Passive) 

92.0 - 95.0 

FIXED SATELLITE (Uplink) 

95.0 -100.0 

AERONAUTICAL MOBILE SATELLITE 
MARITIME MOBILE SATELLITE 
AERONAUTICAL RADIONAVIGATION .SATELLITE 
RADliiNAVXGATtON SATELLITE"" 



1 


A review of Table 3 indicates that most of the frequency spectrum 
above 10 GHz is assigned to the satellite services or the radio astronomy 
services. This does not mean that the FCC and DOC will utilize them, but it 

9* ’l 

does h-ighl ight the potential for use of these frequency bands. In Figure 1 
is shown those frequency segments not assigned for potential use by the 
services listed fn Table2. ! 

I ' 
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4 .0 OTHER PROPAGATION EFFECTS NOT ADDRESSED LN THIS REPORT 


4 .1. IONOSPHERIC EFFECTS 

The Ionosphere generally has a small effect an the propagation of 
radio waves In the 10 to 100 GHz range. Whatever efieects do exist 
(scintillation, absorption, variation in the angle of arrival, delay, 
frequency change and depolarization) arise due to the i ntaraction of the 
radio wave with the free electrons, electron density irregularities and the 
earth's magnetic field. The density of .electrons in the ionosphere varies as 
a function of geomagnetic latitude, diurnal cycle, yearly cycle, and solar 
cycle (among, others). Fortunately, most U.S. ground -station-satellite paths 
pass- through the midlatitude (lowest) electron density region yielding an 
even lesser effect on propagation. Canadian stations may be affected by the 
auroral region electron densities which are normally rhigher. A more complete 
discussion of these effects is included in CCIR Report 263-3 (Ref. 4). 

.A mean vertical one-way attenuation for the ionosphere at 15 GHz for 
the daytime Is typically 0.0002 d8 (Ref. 5), the amplitude scintillations are 
generally not observable (Ref. 6) and the transit time delay increase over 
the free space propagation time delay is of the order of 1 nanosecond (Ref. 
7). Clearly for most systems operating above 10 GHz ilhese numbers are 
sufficiently small that other system error budgets' wliTS be much larger than 
the ionospheric contributions. 
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TROPOSPHERIC DELAYS 


4 .2 


Highly accurate satellite range* range-rate asset position-location 
systems will need to remove the propagation group delay effects introduced by 
the troposphere- Extreassly high switching rate TDMA systems may require 
these corrections in the future. Tne effects arise primarily due to the 
oxygen and water vapor in the lower troposphere. Typical total additional 
propagation delay errors have been measured to be of the order of 8 
nanoseconds' (Ref- 8). 

Estimation techniques, based on the measurement of the surface 
pressure, temperature and relative humidity have been developed (Refs. 8, 9, 
10) which can readily reduce this error to less than 1 nanosecond. In 
addition, algorithms for range (Ref. 11) and range-rate (Rtef. 12) have been 
prepared to reduce tropospheric contributions to satellite tracking errors. 

Since this topic is quite specialized anc generalllly results in an 
additional one-way delay of less than 10 nanoseconds it as net addressed 
further in this report. An overview of this subject and a&iitional 
references are available in CCIR Report 564 Revised (Ref. S3) .. 
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1.0 OVERVIEW 


Earthr-space oath attenuation for systems operating above 10 GHz will 
usually degrade the performance of the path. Consequently, system engineers 
must consider bath the magnitude, frequency and durationof attenuation events 
"along the path?, to realize a workable system. However, because of the random- 
ness of events; in the troposphere only statistical predtrtions can be made of 
path performance. The predictions usually apply on several time scales: 

-long times of the. order -of minutes or more (resulting in an average attenuation) 
and short times, of the order of seconds (termed amplitucb scintillations). 

The scinti Hattons are always present to some degree and arise due to the 
' inhomogeneity of the troposphere. i 

The chief contributors to the average attenuathon along an earth- 
space path are:- ; 

• Gaseous attenuation 

(primarily due to oxygen and W3ter vapor in the 10 to 
1C® GHz frequency range) 

• Rain. 

Some minor" levels of attenuation can be obtained from clouds, fog, 
dust, and sandv The latter contributors have not bseg studied in detail, 
but the results to date will be included to demonstrate the existing level of 
knowledge. Show and ice appear to cause very low leveTerof attenuation, but 
can contribute to wave depolarization. 
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! 

The gaseous specific attenuation, (i.e., attenuation per unit 

i 

length) has been studied and' measured so that its contribution can be 
predicted quite accurately given the path length of the link through the 
medium. - The rain attenuation, however, can only be statistically obtained, 
even though the specific attenuation of rain is now fairTy well knov/n. 

The rain attenuation is usually presented in two formate.. The first format 
consists of the percentage of time the rain attenuation. exceeds a specified 
level (long: tern cumulative statistics) based .on knowledge of- the rainfall 
ir. a giver region, the specific attenuation, and geometric parameters of 
the path. The second format consists of the worst monte statistics (or 
better, the worst: 30-day statistics) based on the local^imatological 
data, the specific attenuation and the geometry of the^aearth- space path. 

Mast- of the rain attenuation -statistics in this report have been de- 
veloped from the common base of the rain rate statistics- Given these rain 
statistics.,, the models develop the total attenuate on staiiisti cs utilizing 
different techniques. Some techniques are easy tO' utilise, while others re- 
quire more, computational time and are expected to be- more accurate. Parts of 
four models which’, have been used in the U.S. have been des:ribed in some detail 
in this report- The- rationale -for utilizing the parti cuter elements of each 
model are described along with the model. However, other models now being 
developed or models not as widely applied are referenced, since they may 
ultimately prove to be more accurate for particular applications. 

Several! approaches to obtaining the rain attenuation statistics 
are presented in; this„ chapter. For any specific design, iit is recommended 
that the simplest approach be applied first and then more sophisticated 
techniques until! the required level of accuracy is obtairat.. The use of 
these approaches - is dependent on the known data. For example, ir no rain 
data is available- for a specific site, utilize the CCIR jGTbbal model. In 
contrast, if long^-term satellite- to-ground attenuation statistics are 
available at another frequency, then only frequency scaling - .will be necessary 
resulting in, a, much more accurate prediction of the link performance. These 
methods a>o described in detail in Paragraph 3 along with specific examples. 

In; the- remainder of this section, data needed for the more detailed 
analyses is also presented. This data forms the basis for the predictive 
model s . 
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2.0 GASEOUS CONTRIBUTIONS* 


In the frequency range from 10 to 100 GHz theavater vapor absorption 
at 22.235 GHz and the oxygen lines extending from 53.5 -±©65.2 GHz are the 
only significant contributors. The absorption lines <are broadened by 
collisions at high pressures (low elevations) thus the specific attenuation 
is altitude deoendent. The vertical one-way attenuation for a moderately 
humid atmosphere (7.5 gm/m 3 } at various station altitudes- above sea level 
is presented irr Figure 2-1. For other elevation angles G (9 >_ 5°) multiply 
the .indicated attenuation by cosecant 0. The range of values C indicated 
in Figure 2-1 refer to the .fine -absorption lines which may be observed between 
high flying aircraft and satellites. (Ref. 2-1). For stairon altitudes not 
shown, a linear extrapolation between curves is recommenced assuming an 
exponentially distributed atmosphere for both water vapor and oxygen. 

The water vapor content is the most variable component of the 
atmosphere. Therefore, for more arid regions a correction should be made 
for the expected values of water vapor. To first order, the vertical one-way 
attenuation at sea-level is given -in Figure 2-2 for 7.5 gn/nf (solid line) 


* The primary references for the material presented in this section are 
Report 234-3, CCIR, Geneva, 1974 and Doc. 5/270-E Workiing Group 5-2, 
19 Sept. 1977. 


• 2-1 


21 


THEORETICAL ONE-WAY ZENITH ATTENUATION (dB) 



1 2 5 10 20 50 .HD- 200 400 


FREQUENCY (GHz) 

FIGURE 2-1 . THEORETICAL VERTICAL ONE-WAY ATTENUATION FR03- SPECIFIED 
HEIGHT TO TOP OF THE ATMOSPHERE FOR A MODERATE HUMID' AMDS PH ERE 

{7.5 G/M 3 AT THE SURFACE) 

A: STARTING' HEIGHTS (KM) 

B: MINIMUM VALUES FOR PATHS STARTING AS INDICATED HEIGHTS (KM) 

C: RANGE OF VALUES FOR THE PATH FROM- THE SURFACE TO $0 KM 
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FIGURE 2-2. THEORETICAL ONE-WAY ATTENUATION FOR VERTICAL PATHS THROUGH 
THE ATMOSPHERE (CALCULATED USING THE UNITED STATES' STANDARD ATMOSPHERE 
FOR JULY AT 45° N LATITUDE). SOLID CURVES ARE FOR A- MODERATE HUMID 
ATMOSPHERE (7.5 G/M 3 AT THE SURFACE); DASHED CURVES FOR VERTICAL 
ATTENUATION REPRESENT THE LIMITS FOR' 0 AND 100% RELATIVE HUMMIDITY 

A: RANGE OF VALUES 

B: RANGE OF VALUES 
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and 0 gm/m 3 (0% dashed Tine) for the US Standard Atmosphere- in July at 45°N 
latitude. Linear interpolation may be used between the carve for 7.5 gm/m 3 
and the 0 gm/m 3 and extrapolated to at least 15 gm/m 3 . T&e el evati.cn angle 
(5° < 6 < 90°} dependence is proportional to cosecant 8. Hhe horizontal curve 
(7.5 gm/m 3 ) in Figure 2-2 can be utilized as a limit for .angles below 10 

elevation angle.. j 

i 

For a station at arbitrary altitude and surface vaxer vapor density, 
a combination of the results in Figures 2-1 and 2-2 is regiired. A suggested 
procedure fallows - : ' Tr • ' 

1. Determine the one-way vertical attenuation the to altitude 
■ from - Figure 2-1 at the frequency of interest 

2. Determine the decrease (increase) in dB ±»elrw- (above) the 7.5 
gm/m?’ curve in Figure 2-1 based on ? linear extrapolation 

-'•between the solid curve and the 0%, 0 gm/m 3 curve. 

3. " Add I)- and 2) in dB and multiply by the cosecant of the 

' ~! ; " r elevation angle 6 (5° <_ '9 <_ SO 0 ). 1 ' ‘ 

. . _ Note: The water vapor content at the surface is usually measured as 

t > ’ r . "« 

relative humidity or the partial vapor pressure. The conversions between 
these units are given in an Appendix at the end of this section. • 

,-f As an example, the GS#C ground station altitude Is 880 meters so 
the vertical one-way attenuation is 0.17 dB at 20 GHz. Using the example 
numbers in the Appendix, namely 30°F, R.H. = 60% and elevation angle * ,47°, 
the- water vapor density is 15 gm/iT? . Extrapolating this'Aalue on Figure 2-2 
yields an additional 0.7 dB attenuation above the 7.5 gm/ifl value. Thus the 
total attenuation is 0.87 dB vertically. Note that this nearly corresponds 
to the 100% RH curve on Figure 2-2, but this is coincidental and indicates 
that .the moisture - level selected for the example is higher than the averages 
displayed in Figure- 2-2. The total one-way attenuation toATS-6 from the 
Rosman ground" station due to gases is computed for a 47-degree elevation 
angle to be 0.87' esc (47) dB = 1.2 dB at 20 GHz. At 30 Glr, the result is 
lower (0.96 dB) for the same conditions since 30GHz is. ne*r the local minimum 
in the gaseous specific attenuation values. 
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3.0 RAIN ATTENUATION 


-3.1 OVERVIEW 

3.1.1 Scattering 

Rain drops both attenuate and scatter microwave energy along an 
earth-space path. From the basic Rayleigh scattering crfteria (the dimen- 
sions of the scatterer are much smaller than the waveTe%th ) and the fact 
that the median rain drop diameter is approximately 1.5 mm, one would expect 
•that Rayleigh scattering theory should be applied in the frequency (wavelength) 
range from IQ GHz (3cm) to 100 GHz (3mm). However, Rayleigh scattering also 
requires that the imaginary component of the refractive fndex be small, 
which is not the case for water drops (Ref. 3-2). Became of this effect 
and the wide distribution of rain drop diameters, the Rayleigh scattering 
theory appears to apply only up to 3 GHz (Ref. 3-2). -Above 3 GHz Mie 
scattering applies and is the primary technique utilized! for specific rain 
attenuation (attenuation per unit length, dB/km) calculations. Mie scattering 
accounts for the deficiencies of Rayleigh scattering andi has proved to be 
the most accurate. 

3.1.2 Drop-Size Distributions 

Several investigators have studied the distribution of rain drop 
sizes as a function of rain rate end type of storm activity. The three most 
commonly used distributions are: 
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fsr. 9 . Laws and Pat sons (LP) 

• Marshall-Palmer (MR) 

s Joss-thunderstorm ( J-T ) and drizzle (J-D) - 

J * 

In general the Laws and Parsons distribution (Ref. 2-3) is favored 
for design purposes because it has been widely tested by caaegarison to 
measurements for both widespread (lower rainrates) and' comsstive rain 
(higher rainrates) at the present time. Also in the higher irainrate 1 R 
regime (R 25 mm/hr) and at frequencies above 10 GHz, the IF values give 
higher specific rain attenuations than the J-T values (Ref*. 3-4). In 
addition it has been observed that the raindrop temperature lis most accurately 
modeled by the 0°C- data rather than 20°C, since for most :h3gh elevation angle 
earth-space links the raindrops are cooler at high altitudes and warm as 
they fall to earth. 

An example- of the measured ;number- distribution of ^raindrops with 
drop diameter as a function of rain rate R (mm/hr) is qi-ven iin Figure 3-1. 

Here the measurements of Laws and Parsons (broken lines, -Ref.. 3-3) and 
Marshall and Palmer (dotted lines. Ref. 3-5) are fitted by exponential 
relation of the form 


where 



mm 


i 



D(mm) 


FIGURE 3-1. RAIN DROP SIZE DISTRIBUTION FUNCTI09 (SOLID STRAIGHT 
LINES): COMPARED WITH RESULTS OF LAWS AND PARSONS (BROKEN 
LINES) AND OTTAWA OBSERVATIONS (DOTTED LINES’) 

REFERENCE: FIGURE 2 OF REF. 3-5 


Note that the urrfts in the equations and figure 3-1 are different. Multiply 
the Np by IQ 5 to convert to the units of Figure 3-1. The number of rain 
drops with diameters between D and D + 6D in' a volume V (cm? } at rainrate 


' ^ = N 0 5D V- 

1 

. As shown in Figure 3-1, the measured data deviates- from the exponen- 
tial relation far diameters below 1.5 mm. flowever, the. larger drops tend to 
dominate the specific attenuation at the higher rainrates of' most concern 
.for the system engineer, and so this deviation tends not to be reflected in 
.the integral over drop diameters utilized in specific attenuation calculations. 

t 

. Joss, et al . , (Ref. 3-6) have found significant variations, of Np 
and ‘A for different types of rainfall based on one year's measurements at 
Locarno, Switzerland. These results are presented inTajslie 3-1; however, tut 
cl-imatic regions where the Joss statistics apply ha ve-noit been determined. 
Therefoie, it appears best to utilize the Laws and Parsonstesults , realizing 
that in certain areas of the U.S. and Canada they have mot been verified. 

VS 


■ ■ TABLE 3-1 

VALUES OF N 0 , A VERSUS RAIN EVENT AS 
DETERMINED BY JOSS, ET AL. (REF. 3-6) 


Rainfall 

Type 

No 

-4 

cm 

, 

i , . 

A'. 

— l 

-am 

drizzle 

0.3 

...05. 2 1 

57R 

widespread 

0.07 

. J ,21 

41R 

thunderstorm 

0,014 

- w - n — 00 ’' 2 1 

30R 
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3.1.3 Specific Attenuation 

It has been determined {Ref. 3-7) that the specific attenuation A 
(dB/km) is related to the rain rate R (mm/hr) by a relation 

A * a(f}R b!f ' b ' 

* 

where the coefficients a and b are functions of frequency.. At this time the 
most thorough calculations of A have been made by Olsen, et al . (Ret, 3-4). 
These calculations extend from 1 to lODOGHz and have been presented in both 
tabular and empirical format for several raindrop distributions and tempera- 
tures. For the U.S. and Canada the 0°C numbers are most applicable (Ref. 3-2) 
Table I of Ref. 2-5 is given below (Table 3-2) for selected frequencies of 
interest in this report. The LP L and LP H refer to Laws and Parsons drop size 
distributions associated with rainrates R from 1.27 tot5E8 mm/hr and 25.4 to 

152.4 mm/hr, respectively. Olsen, st al . have also pro.viided analytic approxi- 
mations for ?.(f) and b(f) which are quite adequate for tee by system engineers 
These are 


a(f) = 4.21 x 10- 5 C f) 2 *" 2 


2-9 < f -< 54 GHz 


- 4.09 x 10’ 2 (f) 0,599 


5T < f < 180 GHz 


and 


, % _0 0779 

b(.f) = 1.41 (f) * 


« 2.63 (f) 


- 0.272 


a, 5 < f < 25 GHz 
25 < f < 164 GHz 


where f is in GHz. 

Thus for 20' GHz 

Pi ~ a(f)R^^^ dB/km 


= 4.21 x 10 " 5 (20) 2 -'* 2 R 1 - 41(20) 0,073 dB/km 


1 117 


= Q.059R * = 2.19 dB/km @ R = 25.4mm/hr. 
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TABLE 



3-2 

R b (dB/km) AS FUNCTIONS g)F FREQUENCY AND 
IN. TEMPERATURE = 0°C " 


■J-D 

— 
- r L . 

' lp h ' 

MP 

d-T 

d-D 

1.14xl0" 2 

1.178 

1.189 

1.150 ( 

1.076 

0.968 

1.41X10** 2 

1.171 

1.167 

1.143 

1.065 

0.977 

1.72xlO" 2 

• 1.162 

1.150 - 

1.136 

1.052 

0.985 

2.82xlO" 2 

: 1.142 

1.119' 

1.118 

1.010 

1.003 

5.30xl0“ 2 

' 1.119 

1.083 

1.097 

0.946 

1.020 

8.61xl0~ 2 

: 1.094 • 

1.029 ' 

1.074 

0.884 

1.033 

0.128 

! 1.061 

0.964 

1.043 

0.823 

1.044 

0.180 

. 1.022 

0.907 

1.007 

0.783 

1.053 

0.241 

0.981 

°1 

0.972 

0.760 

1.058 

0.387 

;■ o:3or 

o'Ms. 

0.905 ' 

C.709 

1.053 

0.558 

: 0.850 

0.794_ 

.0 • 85 1 

„0.602 

035 

d»?ld 

0.809 

6,784 

0.812 

016^1 

LQ09 

0.922 

6.778 

o 

CO 

« 

o 

0.781 

0.674 ( 

0.980 

1.10 

0.756 

0.776 

0.753 

0.663 

0.953 

1.26 

0.742 

0.774 

0.730 

0.637 

0.928 


77 '/ 



_ ' . 1 ^ 119 - 

The value in Table 3-2 for this frequency is C.0626 R =2.34 
dB/km § R = 25.4 mm/hr, an error of 6 %. 

" Based on the Olson, et al , , results (Ref. 3-4), the following spe-= 
cific attenuation- relations are recommended in the 10 to 100 GHz frequency 
range (see Table 3-3). These specific attentuations vriTH be utilized through- 
out this report. j 


TABLE 3-3 

• RECOMMENDED SPECIFIC ATTENUATION APPROXIMATIONS, 


i 

Frequency 

Specific Attenuation , A, (dB/km) 

~ 

Range 

(R in mm/hr,. f in GHz) 

L i 

V ,<^ 


r- I 

v‘ • 

- ‘<y/ 

10 - 25 GHz 

.-.-0 .0779]’ 

A = 4,21 x 10" W- 42 R U;**i(' f ') j 

I. 

L ' *.■=! 

1 

i 

25 - 54 GHz 

A = 4.21 x 10" 5 (vf)) 2 * 42 R .63(f) °* 27i ] 


54 - 100 GHz 

A = 4.09 x 10~ 5 (lf); 0 ' 699 R l 2 - 62 ( f ) * ' J 

j 


The specific attenuations for several of the cimmon earth-space 
bands are shown in Figure 3-2 for raindates from 0.1 to 10 inches/hr (2.54 
to- .254 mm/hr) using the equations in Table 3-3. The 85 and 94 GHz curves 
overlap the 50 GHz data because of inaccuracies in the approximations in 
Table 3-3. More accurate results are obtained from interpolation of 
Table 3-2. 
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SPECIFIC ATTENUATION (dB/Km) 







3.2 


RAIN RATE DATA 


3.2.1 U.S. Sources 

Ire the tLS. the National Weather Service's National Climatological 
Center* prepares, and' maintains extensive precipitation retards obtained from 
Weather Service Offices and over 12,000 observers and agencies. This rain 
data is available in several documents available from the National Climato- 
logical Center.. Several of these publications of interest to the earth-space 
path engineer are:: 

* Hourly Precipitation Data (HPD) 

15 minute rain rate resolution 
published monthly by state 

District of Columbia included in the Virginia HPD 
available about 6 months following date of recording 
$0.4G per copy 

»s Climatological Data (CD) 

1 hour rain rate resolution 

- published monthly by state 

District of Columbia included in the Maryland and 
Delaware CD 

- Washington National Airport WSO is included in the 
Virginia CD 

available about 3 months following data of recording 
$0.40 per copy 

• •’ * « 

9 . Climatological Data-National Summary 

greatest 24 hour rain rate. data 

- published monthly 

available about 4 months following date of recording 
$0.70 per copy 

• Climatological Data-National Summary, Annual: Summary 
one 5 minute rain rate resolution per .month 
available about 18'months following l$st date of 
recording 
$1.80 per copy 

* National Climatic Center, Federal Building, Asheville,. North Carolina 28801 
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• Local Climatological Data (LCD) 

- hourly rain- rate resolution 
published monthly by location 

available about 4 months fo 1 lowing date of" recording 

$0.20 per copy j ’ 

- annual issue also published for each location, 

$0.15 per copy 

• Stomp Data 

published monthly for the U.S. 

- describes type of storm and extent of .damage. 

The Local Cl imato logical Data is available for the 291 stations shown in 
Table 3-4; however,, the Hourly Precipitation Data is available for many 

more stations. 

Examples of the precipitation-related data available in each of 
these publications, are given in Figures 3-3 to 3-7. Comparing the results 
for either the Baltimore Weather Station Office (VISO) at the Airport (AP) 
l or the Belts.vi'lle results, one observes that precipitation- data up to 15- 
minute; resolution is available in the HPD's, while the monthly CD lists 
only the total precipitation per hour (see Figure 3-4). The monthly CD, 
National Summary, Tists the total .precipitation per month in liquid form 
and-the total snaw-or ice pellet depth at most airports. Also included 
(see Figure 3-5) is the nunfcer of thunderstorms recorded during the month. 

In the Annual Summary of the National CD (see Figure 3-6) tha total precipi- 
tation, snowfall ('all frozen precipitation except hailstones) and the amount 
and date(s) of the; highest precipitation rate during the yesr for periods 
of 5 to 180 minutes are given. Unfortunately it only incliriss^one 5 minute 
event per month, f.e.; if there were several significant high; rain rate 
periods per month, only the highest will be indicated in the data. Addi- 
tional techniques to retrieve more data will be described below. 

The Local Climatological Data (LCD) provides the rainfall by hour 
at each of the 291 stations shown in Table 3-4. An example for Asheville, 
NC, is shown in Figure 3-7. In this publication the type of weather is 
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TABLE 3-4 

LOCAL CLIMATOLOGICAL DATA STATIONS 

U.S. DEPARTMENT 01-’ COMMERCE. 

NATIONAL OCEANIC AND ATMOSPHERIC AD4H2SISTRATIOM 

. ENVIRONMENTAL DATA SERVICE _ 

» i 

r 

(Stations £anv.hich Local Climatological Data are issued, asofJanuary I, 1978) 


ALABAMA 

ibC BIRMINGHAM 
■be HUNTSVILLE. 
tbc MOBILE 
«x MONTGOMERY 

ALASKA 

tbC ANCHORAGE 
M ANNETIC 
M>C BARROW 
13C BARTER INLAND 
*bc BETHEL 
«X BLTTLfS 
M* B'C DELTA 
*K COLO BAY 
*bC FAIR RANKS 
AC COLKANA 
kx HOMER 
-6C JUNTAU 
*bC KING5ALM0N 
toe KOLIAK 
■be KOT7ERUC 
*bC MeGHATH 
Jbt NOME 
tU 5T,PAUw ISLA-NO 
UK TALKEETNA- 
SH UNAVAKLEET 
•be VAUbCZ 
■be VAKUTAT 

ARIZONA 

■be FLAGSTAFF 
■be PHOENIX 
■be TUCSON 
■be WINSLOW 
«Sc YOU A 

ARKANSAS 

. «bc FOR’' SMITH 
XX LITTLE POCK 


■ex BAKERSFlELO- 
’tbc BlS**0" 

MK B<-Uv CANYON 
■C EUREKA 
XK FRESNO 
** LONG UCACH 
•ex LOS ANLt UES AIRPORT 
, ac LOS ANGELES Cl VtC CENTER 
•be MT. SHASTA 
Mx OAKLAND 
■be RED BLUFF 
XX SACRAMENTO 
MX SANORLRG 
Sbe SANOICGO 

SAN FRANCISCO 
•be AIRPORT 
K CITY 
•be SANTA MARTA 
■be STOCKTON 

COLORADO 

MX ALAMOSA 
MX COLORADO SPRINGS 
MX OtNVER 
■be GRAND JUNCTION 
MX PUEBLO 




RLXWiDA 

tX lAPAtACHlCOLA 
'DAYTONA BLACK 
M»r TOUT MYTHS 
KE£ JACKSONVILLE 
*t YWtST 
MC AAXELAhO 
•rtbr dllAW 
<RK ORLANDO 
■MC PENSACOLA 
xrx TALLAHASSEE 
Jx TAMPA 
Mat TWEST PALM BEACH 

GEORGIA 

Abe ATHENS 
Abe ATLANTA 
Abe AUGUSTA 
aoc COLUMftUS 
Abe MACON 
*=x ROME 
•he SAVANNAH 

MAW All 

•be HILO 
•ex HONOLULU 
•be KAHULUt 
•be 4JHUE 


CONNECTICUT 

•be BRIDGEPORT 
XX HARTFORD. 

OEL AWARE 

MX WILMINGTON 

DISTRICT OF COLUMBIA 

MX WASHINGTON • NAT7X3NAL AP‘ 
MX WASHINGTON -DtXLSlSJNrLAP 


Mac iBO ISC 
ax aXwtSTON 
■fee tPOCATELLO 

R03MTS 

ac CAIRO 

{CHICAGO 

XX AL (WAY AIRPORT 
•C tOHARE/.IHFORT 

mu: jmounc 

■tac.7£ORiA 
ee .ROCKFORD 
tex? SPRINGFIELD 

IS8DIRKA 

■ac lEVAJiS/ILLE 
■C3C -ORTV.AVNE 
xe awoiANAPOLtS 
me ^OUTMeENO 

tCWJL 

■KX BURLINGTON 
*t* DCS MOINES 
MX -DUBUQUE 
•be SiOOXOTY 
MX -WATERLOO 

■UCKSAS 

Mtx CONCORDIA 
«w tXJOGCClTY 
■OX GOOOLANO 
•be TOPEKA 
•e -WICHITA 

KENTUCKY 

•be LEXINGTON 
•be LOO»S VI UX 

LOUISIANA 

MZK ALEXANDRIA (f) 
Mae BATON ROUGE 
MX .LAKLLHaHLtS 
arc HFw ORLEANS 
<at SHREVEPORT 


■CARiGOU 

■PORTLAND 


RYLAHO 

BALTIMORE 


MASSACIIL^LTTS 
■be BOSTON 

■c Blue hill cas. 

•be WORCESTER 


MIC HI GAN 
MX ALIY HA 

f>LTPOlT 

•be CtTYAIHPCRT 
•be DETROIT METRO AP 
•tx FLINT 
•be CkANfj AAPiOS 
/be HOUGHTON LAKE 
MX LANSING 
K MARQUETTE 
■be MUSKiGON 
•be SAULT SIC MARIE 

MINNESOTA 

■be OULIITH 

• be INTER national falls 
■ tx Minneapolis st Paul 

■oc ROCHESTER 
MX ST CLOUD 

MISSISSIPPI 


MISSOURI 

■be COLUMBIA 
■tx Kansas city 
• be ST JOSEPH 
■be S’. LOUIS 
■be SPRINGF.EI 0 

MONTANA 

•be BILLINGS 
*oe GLASGOW 
■be GREAT FALLS 
■be HAVRE 
■oc HELENA 
•be KALiSPFLL 

MX MLESCTV 

■be MISSOULA 

NEBRASKA 

MX GRANO ISLAND 

•be LINCOLN 

•be NORFOLK 

■be NORTH PLATTE 

■be OMAHA 

■C OITAHA (NORTH) (2) 

mx scottsolufp 
• e VALENTINE 

NEVAOA 

•be ELKO 

■ be ELY 

*tx LAS VCCAS 

■be RENO 

■be W.NNEVUCCA 

NEW HAMPSHIRE 
■be CONCUR U 
■c MT WASHINGTON 

NEW JERSEY* 

ATLANTIC CITY 
■be AIRPORT 

• state marina 

(M NEWARK 
*C TRfNTON 

NEW MEXICO 

■be ALBUQUI RQ’JE 
■bC CLAYTON 
•K ROSWELL 


YORK 

Sac- ALBANY 
uxy H'NGHAMTON 
BUFFALO 
NEW YORK 
■bC* CENTRAL PARK 
■bexJF KENNEDY 1MT-LAP 
.bb LMJUAR01A TlELO 
ROCHESTER 
je- SYRACUSE 

*JORTH CAROLINA 

■•c* ASMEVItLf 
at*. CAPE HATTER AS 
*»C* CHARLOTTE 
•»C GREENSBORO 
MC RALEIGH 
•■bC* WILMINGTON 

rOGTH DAKOTA 
mx eiSV-AACK 
JXX FARGO 
XXXWUJSTON 

oao 

s*C< A*RC**CAIfTCN 
ONOWlATl 
•x A?EE0tfi. 

;«xA!RfORT 
tRte. CLEVELAND 
XXxCXXLAeuS 
•XxOAYTOfl 

~xx VANSFIEUD 
rJ or TOLEDO 
‘KK tTXf^STWN 

OWOM 


OEGCN 

-JbC ASTORIA 
MbtboftNS 
^xx EUGENE 
jxx* COMORO 
w PENCLSTCN 
-*c ?CRtL*NU 
XX SALEM 
ioc< SEXTCft SU>*MT 

'RafjciSLWts 

eh. GUAM 
*X*JOWON 
<*c KOROR 

KWAJALEIN 
■XX UVUUHO 
•Mi P*oOPAGO 
‘bOJFO Wt 
awe tRLK,moCnj 
• tCxVIAKE 
•■■4XYAP 

*Y»<SYLVAN1A 

Jate AUJEN7 OMi 

•sex AVQCA.WILKESEARRE SCRANTON AP 
^txERiE 

i *cc HARRISBURG 
^JOC HHLAC*LP «A 
PITTSBURGH 
**tC AI fCORL 
tc OTY t 

1-CcWIEUA* SPORT 

FOEISIWD 

-X BlOCKISLAAD 
■XT PROVIDENCE 

•SOTH CAROLINA 

CHARLESTON 
£4£e AtRtCRT 
■■ CITY 
■^‘ONJWBiA 

U3t£. GRCENYUll SPAATANOUHG 


SOUTH DAKOTA 

•be AIBRIFCN 
4t* HURW 
■be RAPiticiTv 
■ex Si (Art FALLS 

UrJASSCE 

«■ hRSTOL 
•«. ciwnw»yGs 
■w ICJOxviilC 
UK MJ.R»S 
•tc NAShvillT 

•c CAKKiOCJ 

nxAS I 

MX ABILFNf- 
•X AMABILLO 
KK AU*TIN 
■be BROl^Vll.r 
*c comuson st 

M< (VILLAS I Lt« ; . AVI 
rCLRIC 
•x (LF'60 
■c CSUftS’Y. 

•he lOtRlON 
*rx LUxYav 
•X v*»:XA.<T 
•< KVITAIIHUR 
xe SananAlO 
ac S/NAfRWO 

■be VI CTO ‘'I A 

MK WACO 

■be .1 CM TA TAILS 

UTAH Y 


4b( BURJ'GTON 

VI ROMA 

■be LtrO*A«G 
•X NORFOLK 
xx RiCre t>jj 
XX ROAiNCXJ. 

*u V.ALlCTS ISLAND 

WASMNGIOh 

•ec Ot.YJ.rtA 
XT QLIUAYUTE AIRPORT 
•be SCAtrietA O* AAI* 
■C SEATTLE U - '. AN SI T t 

xk sro*uci£ 

XX STA\M.rePAjS 

i- WALLA V\ Tu_.\ 

XK YAVJfM 

VESTUOES 

XX SA.NXJATL P R 

WEST VIRG-NI A 

•X WCKLfY 
■re CkahlESTln 
■be Elkins 
jx •ajntji*.ttn 
•». PARVJ-RSPU*IG 

IMSCCYCJN 

XX GRIYNPAY 
JX LALRCAM. 

Jtv tVUAbCIL 
■be MLW.UMLt 

WYOKtNG 

•X CAVER 
•DC Cltvuw 
XX LAND* (I 
■tX Si 4 RICAN 


L Monthly umnrr UKlurtn rml»N» i»iivm,l*em 


t AtYmM Su-rnwy 


AH 


httxn pbOftKY LkjI <-nmtM<IO<)l» 4 l (NtJ I’ 'j\ PM ytJI *rxxul «l It *> 
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HOURLY PRECIPITATION 
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FIGURE 3-5. AN EXAMPLE OF NATIONAL SUMMARY 
OF CLIMATOLOGICAL DATA ISSUED MONTHLY 
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provided so that one can ascertain if the rainfall is from; a thunderstorm or 
a general wi de-coverage weather system. The water equivalent of the snow is 
shown in the hourTy precipitation data. An example of tise LCD rainfall data 
for Asheville, NO, for the month of August 1975, is shown fn Figure 3-7. 

Note that the same information is available' on the Hourly Precipitation Data 
records but that the type of rainfall event is not noted iin - the latter. 

Finally' the National Climatic Center prepares a '"Storm Summary on a 
monthly basis. This information is of little value to system engineers since 
it emphasizes the damage done by the storm rather than "the physics of the 
storm. For example, the most severe rain event in Asheville, NC, in 1975 
occurred on August 24; however, this is not indicated in Hie Storm Summary 
because it apparently caused no significant damage. 

If more information is desired regarding higher stain rates associated 
with .finder storms it can be obtained for most first-order Weather Service 
.Office (defined as those offices manned by Weather Service personnel } sites. 

. i 

These sites generally have both tipping bucket and universal weighing gauge 
precipitation monitors. The tipping bucket gauges generally accumulate the 
number .of 0.01 inch, precipitation events in a day which is utilized to 
collaborate with the accumulation in the other gauges. Hcvever, some tipping 
bucket gauges employ a readout strip chart (triple register chart of opera- 
tions recorder register) similarto that shown in Figure 3-8. By estimating 
the time between tips the rain rate may be estimated. The location of 
those stations having triple register charts was not available from the 
National Climatic. Center. 

The universal weighing gauge is also capable of providing rain rate 
information and is the main instrument utilized to providf. the 5-minute to 
1 hour precipitation data. This measurement is accomplished by reading 
directly from the 24-hour strip chart on the gauge. An -scamp! e of one of 
these strip charts is shown in Figure 3-9. These charts are available 
dating b4ck a'boiLt 10 years from the National Climatic Center for 25 cents 
per chart. By measuring the slope of the line, the rainnate to at least 
5 minute resolution may be obtained and even 1-minute -rain rates may be 
inferred from some charts. It appears that these charts are the best 
source of information for short duration rate ‘data. 
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WS FORM 102tiC U.S. DCPAHTMCXr OFJ2 d«KCR«-^OAA WAnSMAL 

12 INCH DUAL OR 6 INCH SINGLE TRAVERS? - 24 HOUR si*nd*rd 

D*yl«,h» ( } 


RECORD OF PRECIPITATION 


MONTH DAY YEAH 


OBSERVER __ 


snvnor* 


ER ^ CHART ON ?'3*% #06r 2f 7-T 

ro FF 8*9.% & 75* 



FIGURE 3-9. AN EXAMPLE-OF -A UNIVERSAL WEIGHING GAUGE STRIP CHART 
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The last automated rain gauge utilized by the U.S. Weather Service is 
F the Fischer-Porter gauge. This unit is a weighing gauge whicn punches a paper 
strip chart in a BCD format every 15 minutes. The gauge may be set to record 
every 5 minutes, but that resolution is generally not utilized by the Weather 
Service. The gauge records to only the nearest 0.1 inch. 

An example of how intense rain rates may be measured i*s now given. 

Note that from Figure 3-6 the most intense rain rates (0.38 inches in 5 minutes) 
in the Asheville, NC, WSO occurred on August 24, ending at 1658 Eastern Standard 
Time. This occurred during a thunderstorm (see Figure 3-7} but it was not the 
most rain in a 24-hour period, which occurred on August 17. The amount of 
precipitation between 1500 and 1700 EST on August 17 is noted; in the LCD in 
Figure 3-7. However,, the most data appears directly in the gauge readout shown 
in Figure 3-9. Clearly the rainrate just before 4 PM was -more than 4.56 inch/hr 
(116 mm/hr) for the first several minutes. Interpolation .yields a rate approach- 
ing 150 mm/hr. for 2 minutes. Another example of a cloud .burst is shown in 
Figure 3-10. Herein rain rates approach 300 mm/hr (12 inches^hr) at 8 PM and 
contributed to the airline crash at this airport at this tins- Clearly the 
fcattenuation at a ground station would be significant for this- ^ype of 2 mi-nute 
p.vent (0.00038% of a year). 

• ■'•v. 

3.2.2 Canadian Sources ~ 

The Atmospheric Environment Office* prepares several documents contain- 
ing rain and snow precipitation data. These documents are. 

e Monthly Record-Western Canada, Part 1, 

Provinces of: 

British Columbia 
Alberta 
Saskatchewan 
Manitoba 

$21.60 foreign per year 
$ 2.00 foreign per issue 


* Head Office, 4905 Dufferin Street-, Downsview, Ontario M3H 5T4, Canada 

t Available from: Supply and Services Canada, Publishing Centre, Hull , Quebec, 
j. KDA0S9, Canada, Make checks payable to Receiver General for Canada. 
Canadians should request domestic price schedule. 
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WS FORM 1028C us mwimiwT or Ctwmec-JIMA NATiwm wrniinijBuaL^ 
12 >NCH PUAL OR 6 INCH SINGLE TRAVERSE - 24 HOUR 1 

RECORD OF PRECIPITATION 1 
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FIGURE 3-10. AN EXAMPLE OF AN INTENSE TMi! EVENT 
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Monthly Record-Northern Canada, Part 2 
Territories of: 

Yukon 

' Northwest 

$12.00 foreign per year 
$ 1.00 foreign per issue 

» Monthly Record-Eastern Canada, Part 3 
Provinces of: 

Ontario 

Quebec 

Nova Scotia 

New Brunswick 

$21.60 foreign per year 

$ 2.00 foreign per issue 

e Canadian Weather Review 

published monthly 

covers about 250 surface stations throughout Canada 
$5’. 40 foreign par year 
$Q.50 foreign per issue 

available about one month following the dada- of recording 

e Snow Cover Data for Canada 

- published once per winter 
covers all provinces/ territories 

- - $1.00 foreign per issue 

- available June/July of each year 

The data in the Monthly Reviews (available about four months following 
recording) is of most importance to the earth-space path engineer. As shown in 
Figure 3-11, the rainfall, snowfall and total precipitation 3H& given for aach 
day of the month. The Monthly Summary table indicates the nuiiber of thunder- 
storms, etc., and the recording rain gauge data for selected cities is given. 
These are the maximum amounts for the duration periods indicated on the date of 
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^occurrence. In addition, the number of hourly periods with rainfall accumu- 
lations between GL01-Q.09, 0.1-0.19, etc., inches is noted, obese data are 
obtained from tipping bucket rain gauges measuring in increments of 0:01 inches. 

The tipping bucket rain gauge data is available for many more 
Canadian locations. The charts from these gauges are available upon request 
from the Climatological Recording Services Branch of the Hea# Office in 
Downsview, Ontario, at a nominal charge. 

3.2.3 Worldwide Sources 

Many -countries prepare meteorological data similar ijp the US and 
Canada. Many of Idles e are on file at the National Weather Service Library, 

Room 816, Gramax Bldg., 13th Street, Silver Spring, MD. .One document, the 
Monthly Climatic Data for the World does list the number of days per month a 
station receives more than 1 mm of rain and the total rainfall per month. 

The data is coarse and can only provide a general indication of the precipi- 
tation climate. An example is shown in figure 3-12. , T his document is available 
for $0.70 per monthly copy about 9 months following the recording date from 
the National Climatic Center. 
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3.3 EFFECTIVE PATH LENGTH 
3.3..1 Overview 

l 

Given the specific attenuation A at a given frequency and rain rate, 
the total attenuation along the earth-space path may be calculated given the 
length of the path: through the attenuating medium. The effective path length 
provides the coupling between the specific- attenuation and the’ total attenua- 
tion. However, because the attenuation varies with distance along the path an 
integral must be calculated. The values of the point rain rate R along the 
path at any time t are not known so the concept of the effective path length 
L g has been introduced. At any given instant t it is defined as 

L e (f,s) A tf,R(s=0) ] . At [f.sj 

f Sl 

- jA[f,R(s)]ds 
s — 0 

where L (f,s) is the- effective' path length at frequency f, time t, and along the 
earth-space path s is zero at the ground station and R(ssfl’) is the point 
rain rate at the ground station. A-j- is the total path attenuation along the 
path s at time t, and s t is a .point along the path above which A is negligibly 
small. Implicitly Sjis a function of the elevation angle 9 and the 0°C isotherm 
height where ;the precipitation exists as ice, and R is a function of s and t 
because the rain rate varies along the path s as dense rain cells pass through 
the link. 

The value oif the integral has been measured by direct observation of 
satellite beacons (Ref. 3-8) and radar diagnostics of rain systems (Ref. 3-9). 

3.3.2 Elevation. Angle Dependence 

The path length s of the attenuating medium may 'be closely approxi- 
mated by the cosecant 6 relation. This assumes a pi anegeonstry which is nearly 
the case for typical' heights to the 0°C isotherm (less than 5 km, see 
Figure 3-16) and elevation angles above 10 degrees. 
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3.4. RAINFALL MODELS FOR THE U.S. AND CANADA 

The information in Sections 3.1 through 3.3 will allow one to compute 
the attenuation along an earth-space path given the point rain rate R. In this 
section four modeTs are presented to provide the cumulative statistics for R. 
The models described here are: 

- 1 G" 

a) the Global Rain Model (Ref. 3-12) 

b) the Rice-Holmberg Model (Ref. 3-13) 

c) the Dutton-Dougherty Model (Ref. 3-14, 3-15) 

d) the Lin Model (Ref. 3-16). 

The results presented herein are only those portions of the model which apply 
to calculation of R and its distribution along the path. In Section 3.5 the 
total attenuation is calculated. 

Other rain and rain attenuation prediction and scaling models have 
been published and are in the process of being developed. These include: 

e) the Synthetic Storm Model developed at the Virginia Polytechnic 
Institute and State University, R. F-. Persinger, W. L. Stutzman, 
R l - E.. Castle and C. W. Bostian, preprint;; and R. R. Persinger 
and W. 1. Stutzman, VPI & SU Interim Report 1S78-1, June, 1978. 

f) the Frequency Scaling Model developed at COMSAT Laboratories, 

F. Tseng and G. Hyde, EASCON 78, Arlington* VA, Sept. 25-27, 

1578, pg. 396. I ' 

g) the Storm Cell Model developed at Ohio State' University, 

L. R. Zintsmaster, ElectroScience Lab. States Report 2374-9, 

M*y 1972. 

h) the Location-Dependent Attenuation Model developed at the Applied 
Physics Laboratory, J. Goldhirsh, Proc. UR5I, Comm. F, La Baule, 
France, 28 April - 6 May 1978, pg. 305; anff Radio Science , 

vol. 12 , Sept - Oct 1977, pg. 741. 

These four modeTs should be consulted, as appropriate, where they appear to 
extend the results of the more frequently utilized models.. 
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3.4.1. 



—This modeT ts currently being proposed zo supplement and update the 
current CCIR model described in Report 563. The climatic regions are defined 
and then the distribution curve for tne surface point rain rate data in each 
region is developed'.. The five rain rate regions for the tf.S., and Canada are 
shown in Figure 3-14. Therein rate distributions for each, of the regions 
are show in Figure 3-15. The adjacent geographic distributions are used as the 
upper and lower bounds for the point-to-point and year-to-year variation of the 
observed annual distribution function. A table of values corresponding to 


Figure 3-15 is shown in Table 3-5. 



is utilized tG account for the path averaged rain rate compared to the 
ground station rainr rate, where D - H/tan-3 is the horizontal! projection (basal) 
length of the path- to the height H of the 0°C isotherm (see Figure 3-16) and 
i (D) and <5 (D) shown in Figure 3-17 are the parameters for the- path averag- 
1 ing model. The path averaged rain rate is then 


for a given probability of occurrence. The effective path length is - H csuj 
For ground stations above one kilometer, it appears desirable to subtract the 
ground station altitude from the value of H given in Figure X- 15. and use this 
value for the effective path length. . _ ^ 


• 3-29A 


52 






FIGURE 3.14A. GLOBAL RAIN RATE CLIMATE REGIONS 
(FROM REF. 3-12) 
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RAIN RATE DISTRIBUTION VALUES (mm/hr) 
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3.4.2 The Rice-Holmberg Model (Ref. 3-13) 

The Rice-Holmberg model is the basis for the CCIR rainfall model used 
in Report 563, In the R-H model the rain events are divided into two groups. 

The first group applies to convective storms (thunderstorms). The second group 
includes all other rain events. The model yields the amount of time T t (R). that 
a t - minute average rainfall rate exceeds a given value R. 

For one-minute time intervals (the motivation for selecting a one-minute 
average is based on weather data and statistics (see Figure 3-20)). 

T 1 (R) = M {0.03fexp(-0.03R) + 0. 2(1-0) [exp(-0.258R) 

+ 1.86exp(-1.6'3R)]} hours 

.where M is the total annual rainfall in millimeters, and 8 is the ratio of the 
rainfall during thunderstorms to the total rainfall for the year. To utilize 
this formula, M the average annual rainfall, must be determined (not including 
snowfall) from the weather records for a given location. The Annual Summary of 
Climatological Data (Figure 3-6) indicates, for example, that the normal pre- 
cipitation for Baltimore averaged from 1941 to 1970 is 1028 mm minus the snowfall 
The snowfall data is obtained from reviewing many years of the data and assuming 
that 10 mm of snow equals 1 mm of liquid water. 

The ratio of the thunderstorm rainfall to the total rainfall, 8, is 
obtainable from the local Climatological Data for the location. For example, 
referring to Figure 3-7, in August 1975 in Asheville, NC, 2.97 inches of rain 
was associated with the 10 thunderstorms during the month. The total precipita- 
tion for August was 3.63 inches, so 8 = 2.97/3.63 = 0.82. When averaged over the 
entire year this number would be much lower. 

Waps showing M and 8 are shown in Figures 3-18 and 3-19. These numbers 
may also be used in the one-minute average rainfall formula given above to obtain 
the number of hours per year during which the rainfall rate exceeds R. The per- 
centage per year is obtained by dividing by 87.66 (there are 8766 hours per year) 
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FIGURE 3-18. MEAN ANNUAL PRECIPITATION IN INCHES IN 
U.S. AND CANADA (1 inch = 25.4 mm) ■ , 
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FIGURE 3-19. THE PARAMETER 0 IN THE RICE-HOLMBERG MODEL 
OVER THE U.S. AND CANADA 
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FIGURE 3-1 8A. MEAN ANNUAL PRECIPITATION IN INCHES 
{V inch = 25.4 mm) . 



FIGURE 3-19A. 


THE PARAMETER 8 


IN THE RICE-HOLMBERG MODEL 
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3.4.3. The Dutton-Dougherty Model (Refs. 3-14, 3-15 ) 

This model is also referred to as the modified Rice-Holmberg model. 

The model yields an improved estimate of T^(R) , i.e., the number of hours, that 
the t - minute averaged rain rate will exceed R at the surface. Usually t - 1 
minute is considered to be. representative of the instantaneous, rain rate. 

In the Dutton-Dougherty model the percent of an average year P^R) that 
t - minute rainfall rates R are expected are grouped into three rainfall rate 
regimes, namely 

P t (R) = 0.0114 T lt exp(-R/R lt ) R > 30 mm/hr 

P (R) = 0.0114 T st exp(-/R/R s t) 5 < R < 30 mm/hr 

P t (R) = 0.0114 (T lt + T t ) exp (-R/Rp R < 5 mm/hr 

Here the R^, R^, and T 2t are linear combinations of M, & and D = 24 + 3M hours 
(Ref. 3-13) by the relation 


R 

T 

R 


it 

2t 

I 

t 


= a,*M + a 0 ^B + a 3t D + a t+ ± S 


it 


2t 


4t 


= b U M + b 2t 4 s 2 


b 3t M t b 4t e + b 5t D + b 6t ± S3 


where S ls S 2 and S 3 are the standard errors which can be shown to be small 
compared to R lt , T 2t and Rj. (Ref. 3-17). The Natl. Tel ecom. Inform. Ag. will 
make available a well annotated computer program (DEGP76) to do these calcula- 
tions (Ref. 3-18). ‘This program does not include all of the smoothing algori- 
thms currently being added to an updated program. Jhe parameters T lt - M/R ]t , 

T and T . are found from the boundary condition equations between the regimes 
st it 1 

at 5 and 30 mm/hr. 

Dutton (Ref. 3-15) has extended the model to include the variance of R 
for a given T t (R). This allows calculation of confidence limits on the average 
cumulative statistics. These calculations are not included in the DEGP76 program 
and require additional regression analyses of data from several locations within 
the same climatic zone. In addition, Dutton (Ref. 3-15) has demonstrated how to 
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obtain a zona-lly averaged rain rate distribution for a climatic zone (similar 
to the curves in Figure 3-15) and it is suggested that these analyses be utilized 
to define the zone boundaries. 


3.4.4 The Lin Model (Ref. 3-16 ) 

Lin has observed that the rain rate probability is approximately 
lognormal. Namely, 


P(R)^P o (0) < Js erfc 


In R-j - In P m 
____ . 


where ’P(R) is the probability that the surface point rain rate exceeds rain 

rate R, P q ( 0) is the probability that it is raining, is the instantaneous 

rain rate, R is the mean value of R. and S D is the standard deviation of In R. 
m i K 

The value of P (0) is usually calculated for a rain rate of 0.25 mm/hr 
o 

or more. This selection is made because 0.25 mm/hr has little effect on links 
below 60 GHz and most rain gauges can not measure more accurately than -0.25 
mm/hr. In addition, since the Weather Service measures rain accumulations 
during a period of one hour, the effect of integration time has been experimen- 
tally estimated to be 

P(R >_ 0.25 mm/hr j T £ 1 min) % %P(R >_ 0.25 mm/hr J T = 1 hr) . 

This relation allows one to estimate P (0) from the HPD's shown in Figure 3-4, 
the LCD's shown in Figure 3-8 or the Canadian Recording Rain Gauge Data as 
shown in Figure 3-13. This formula is dominated by the low rain rates and should 
not be applied to higher rates where convective storms lasting several minutes 
yield’ lognormal statistics. 


The initial value selected for R m should be obtained from the total 
rainfall for the year (precipitation minus snowfall) divided by the number of 
hours the rainfall exceeds 0.25 mm/hr (0.01 inch/hr). This data is found in the 
HPD's, the LCD's or the Canadian Monthly Review. The final values for R m and S R 
are obtained by a least-squares approximation that minimizes the differences 
between the data points and the lognormal approximation. Typical values for 
Merrimack Valley, MA (Climate Region D, see Figure 3-14)"are (Ref. 3-16) R m ~ 
1.23 mm/hr, S R a 1.34, P Q (0) « 0.033 = 3.3%. 
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The instantaneous rain rate R-j is assumed to be averaged over a 2- 

f 

second interval which corresponds to the volume in the' main Fresnel zone 
divided by the raindrop terminal velocity. Measurements of the rain rate 
distribution dependence on the gauge integration time have been made in the 
range 1.5 sec <_ t <_ 2 min. for at least one location. The data is shown in 
Figure 3-20 from Ref. 3-19. 


The probability P o (0) is the probability of rainfall at a point. 
However, intuitively one would expect that for a path of elevation angle e 
that rain may be occurring along the path but not at the gauge. This has 
been compensated for in the Global model by the factor r. In the Lin model , 
the P 0 is modified to be P 0 (D), where D is the horizontal projection of the 
path. Experimental data in Florida has yielded an empirical relation- (Ref. 
3-20) 


P 0 (D) ~ 1 


1 - P 0 (°) 

1 0.014 



0.014 

+ 21.5 J 



where D is in kilometers. Thus for earth-space paths, the P 0 (P) value shoal-! 
be used instead of P o (0). 

Later models developed by Lin have recently been published and are 
also recommended for use here. Measurements at US locations have shown that 
the total attenuation = (AL)/(1+L/U(R)) , where A is the specific attenua- 
tion, L is path length through the storm, R is the 5-minute point rain rate 
in mm/hr and ]T(R)=(2636/(R-6.2))km (see Refs. 3-21 and 3-22). Typically Lin 
finds that storm height are 4 km, so L-(H-G)csc 9, where G is the ground 
station altitude. 


3.4.5 Worst-Month Statistics 

Designers of communications systems who must meet "any month" per- 
formance criteria v/ill need a model for the worst-month or worst 30-day sta- 
tistics. These statistics are the average percentage of time within the worst 
month which comes once per year when a threshold is exceeded. 

Recently Crane and DeBrunner (Ref. 3-23) have applied an exponential 
model to fit the tail of the distribution curve. Using the notation in Ref. 
3-23, the probability of exceeding the worst-month threshold value is P(xhj), 
where h denotes the worst month and j denotes the threshold rain rate value 
related to the path attenuation margin budget and other system parameters. 
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FIGURE 3-20. THE INSTANTANEOUS RAIN RATE DISTRIBUTION 
VERSUS GAUGE INTEGRATION TIME 
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The distribution for N years is 

P( x hj) 12 N ~ ^oj exp(-X[ ) j/C 1 j ) 

where x^j is the probability of exceeding j during the worst month. Invent- 
ing the equation yields 

x hj = ^i0 ] n (12 N C 0 j) = C lt ln(C 0 j/P(x[ 1 j ) ) . 

An example of this distribution is shown in Figure 3-21. The values C 0 j ana 
C X j may'be obtained from the figure, e.g., C 0 j = 0.2 and C r j = 7.5-x 10" 5 . v 

For the case where ;i = 1 year, a bound on C 0 j may be obtained from 
12 1 n ( 12 C o1 ) 

fl c o 3 - ° J " -l 4 - 41 

where M (M > 3) is the number of months during the year that rain typically 
occurs. If M < 3, the use of the exponential distribution model should be 
questioned. 

The data in Figure 3-21 is best obtained from the CD (Figure 3-4) 
where the number of days with heavy precipitation is recorded or from Table 7 
of the Canadian Monthly Record (Figure 3-11). Ultimately, the rain gauge • 
charts for the dates and locations identified above, will be required to ob- 
tain the high rain rate data. 


6? 
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3.5 TOTAL RAIN ATTENUATION PREDICTIONS 

•3.5.1 General Approaches 

An estimate of the rain attenuation cumulative statistics may be 
determined in several ways. The optimum way depends both on the amount of 
rai.n and/or attenuation data available and on the level of sophistication 
desired. In general the cumulative statistics can be calculated over a 
period of a year with 99% error bounds, so that only 1 year in 100 will be 
the exception. The system engineer, based on the path margin available 
for attenuation, can then judge how the attenuation will affect the per- 
formance and availability of the system. Depending on the data base available,, 
the designer may compute the cumulative statistics in several ways; however, 
it is recommended that the simplest calculations be carried out firsu to 
give a general feeling for the statistics and also to act as a check on the 
results of more sophisticated calculations. 

The several flow charts in Figures 3-22 to 3-24 will assist with 
deciding which calculations are to be pursued. The steps are numbered to 
allow easy reference with the discussions to follow. 

Generally the yearly cumulative statistics are desired. The worst 
30-day or worst-month statistics were discussed in an earlier section. Again 
it is emphasized that only the rain statistics are desired, since snow has 
a much lower specific attenuation. The examples used here are for synchronous 
or near-synchronous satellites with relatively stable orbital positions. 

3.5.2 Analytic Estimates of Attenuation Statistics (Figure 3-22) 

Referring to Figure 3-22, the required given parameters are the 
ground station latitude, longitude and elevation; the satellite orbital 
position which allows calculation of the elevation angle; and the operating 
frequency. Step 1 is to determine the rain zone climate region from Figure 
3-14 and utilize the appropriate one-minute average 'surface rain rate curve 
shown in Figure 3-15. Given the latitude, elevation and season of the year, 
the averaged rain rate along the link is computed using Figures 3-16 and 3-17. 
In Step 2 the effective path length is computed from L g = H cscG and the total 
attenuation is then computed from Table 3-3 and L e . 
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FIGURE 3-22. ANALYTIC ESTIMATE OF CUMULATIVE 
RATE AND ATTENUATION STATISTICS 


RAIN 
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STEP 11 



SPECIFIC ATTENUATION 


FROM aR b 

STEP 13 

COMPUTE: 



RAIN RATE 


EFFECTIVE PATH LENGTH 



COMPUTE: 

CUMULATIVE 
ATTENUATION STATISTICS 



FIGURE 3-23. ATTENUATION STATISTICS GIVEN RAIN RATE STATISTICS 
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STEP 17 




FIGURE 3-24. 


CUMULATIVE ATTENUATION STATISTICS GIVEN LIMITFD RAIN RATE AND ATTENUATION STATISTICS 














3. 5. 2.1 Example: Analytic Estimate Procedure (Figure 3-22) 

3.5.2. 1.1 Global Model. For the case of the 20 GHz ATS-6 link from Rosman, NC, 
latitude = 35 degrees N, longitude = 82 degrees W, elevation = 880 meters, the 
elevation angle = 47 degrees, so it is located in rain region D. The 0°C iso-^ 
therm height H = 3.7 km so the horizontal path length D - H/tan 9 = 3.7/tan(47 ) 

= 3.45 km a 3 km. The rain rate average parameter r = y(P) R “ 1*6 R > 

so the averaged rain rate = r R - 1.6 R 0 *® 7 . This factor is used to, modify the 
ordinate of Figure 3-15, i.e., for a surface rain rate of 50 mm/hr, the averaged 
rain rate is 48 mm/hr. The effective path length is L g - H csce = 5.1 km. The 
specific attenuation (from Table 3-3 ) is 

A =0.06 R ave 1,12 dB/km 

where the path averaged rain rate R ave is utilized. The total attenuation due .to 
rain will exceed 

AL e = 0.06 (48) 1 - 12 (5.1) dB = 23.4 dB 

for 0.01% of a year (53 minutes) using the Global Model. The lower bound for 0.01% 
of the year is computed from curve E of Figure 3-15 and the upper bouiid from curve 
B of Figure 3-15. 

3. 5.2. 1.2 Rice-Holmberg and Lin Models. If a more accurate estimation is desired, the 
Rice-Holmberg model may be applied (see Step 5). The M and B paramters may be computed 
by averaging several years records from a nearby Weather Station Office or from the 
maps in Figures 3-18 and 3-19. Using the data in these Figures for the Asheville, 
NC-area, M 1150 mm and g = 0.3. For the surface rain rate of 50 mm/hr 

T ! (50) = 1150 | (0.03) (0.3) exp [(0.03) (50)] 

+ (0.2) (0.7) exp [(-0.258) (50)] 

+ 1.86exp [(-1.63) (50)] } 

=2.3 hours = 0.026% of the year. 

Clearly^ this is higher than the 0.01% predicted by the Global Model. 

The Lin model would use the 5-minute rain rate data from the Weather 
Station office. This has not been done in this report, however, data for Palmetto, 

GA (Ref. 3-22) yields T^50) = 1.1 hours = 0.012% of the year. 


*Some compensation for the ground station height has been considered by 
rounding 3.45 km down to 3 km. 
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3.5.2. 1.3 Dutton-Dougherty Model. Because of this disparity it appears 
worthwhile to check these estimates with other models (e.g., the Dutton- 
Dougherty model, Step 8). The results for all three models along with 
the measured ATS-6 data are shown in Figure 3-25. The Dutton-Dougherty 
model best approximates the measured 20 GHz data from ATS-6. Unfortunately, 
the Dutton-Dougherty model also includes the clear attenuation of approxi- 
mately 2.5 dB. This value appears high (see calculations in Section 2 
yielding 1.2 dB) and so the Dutton-Dougherty model results should be shifted 
by about 1 dB. 

3. 5. 2. 2 Discussion of Example . In this case the more difficult model to 
apply yields the most accurate results. The Rice-Holmberg model fits quite 
well for low values of attenuation but begins to deviate at values of 
probability in excess of 0.02% or less. For percentages near 0.1% the model 
tends to underestimate the data. The Global model also tends to under- 
estimate the statistics of the attenuation. Similar estimates made for 
White Sands, NM (Region F) also indicate that the attenuation is under- 
estimated with the Global model. 

3.5.3 Attenuation Estimates Given Rain Rate Statistics. (Figure 3-23) • 

If the rainfall is measured by a weather service for a period of 
10 years or more near a ground station site, the rain rate statistics may 
be generated from the LCD's or HPD's supplemented by the universal weigh- 
ing gauge charts (see Figures 3-9 and 3-10) or a tipping bucket chart as 
shown in Figure 3-8 during the heavy rain events. It is desirable to 
obtain the data in one-minute rain rate averaged increments. This is close 
to the instantaneous rain rate needed for the statistics (Ref. 3-19). 

The rain rate statistics are now averaged along the link path 
using the relation R avg = r R (Step 12A), or the point rain rate may be 
utilized directly in the specific attenuation formula (Step 12). If the 
averaged'- rain rate is utilized, then the effective path length is L g = H csce. 
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PERCENTAGE OF YEAR ATTENUATION EXCEEDS ORDINATE' 


POINT RAIN RATE (mm/h'r) 


20 GHz, ATS-6, ROSMAN, NC 
A GLOBAL MODEL 
O RICE-HOLMBERG MODEL 
0 DUTTON-DOUGHERTY MODEL 


□ ATS-G DATA OBTAINED 
BY DISTRIBUTION 
EXTENSION TECHNIQUE 
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ATTENUATION (dB) AT 20 GHz 


FIGURE 3-25. CUMULATIVE ATTENUATION PREDICTIONS (3 MODELS) 
AMD’ DATA FOR 20 GHz SIGNALS AT ROSMAN, N.C. 


E 









If the point rain rate is utilized in the specific attenuation calculation, the 
effective path length is found in Figure 3-13 or the formula L g (R) = 12.82 R c ‘ ' 
should be utilized. Since the effective path lengths in Figure 3-13 have been 
measured, it appears more accurate to use the experimental relations. The 
results of these calculations will be the total attenuation at the operating 
frequency. These attenuations are then plotted on the parallel to the rain 
rate axis similar to Figure 3-25. 

3. 5. 3.1 Example of Attenuation Statistics Given Rain- Rate Statistics. 

(Figure 3-23) . At Rosman, NC, the rain rate from the near bucket 

was measured over a period of six months (July to December, 1974). The 
resulting distribution is shown in Figure 3-26. The attenuation values for 
a given probability of occurrence are then computed using the specific 
attenuation and the effective path length relation L g = 12.82 R 0 ’ 3 . The L g 
term causes the attenuation estimate to deviate from the rain rate distribu- 
tion for high values of rain rate. 

3.5.4 Attenuation Estimates Given Limited Rain Rate and Attenuat ion 

Statistics (Figure 3-24) 

3. 5. 4.1 Distribution Extension . Referring to Figure 3-24, the amount of given 
data is more extensive than the other methods mentioned, except that both the- 
rein rate and attenuation statistics were measured over a limited period of 
time (typically less than a year). It is desirable to extend these statis- 
tics to a one-year period through a procedure called distribution extension. 

Consider the rain rate distribution first. Two situation are likely 
to occur. First, if the measurements were made for a total time T only 
occasionally during a time period T x to T 2 , and one desires to extrapolate the 
measured distribution to the total time T 2 - T 1 . The second situation occurs 
when the measurements made between Tiand T 2 are to be extrapolated to a larger . 
time period T 0 to T 3 which includes Ti and T 2 . 
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PERCENTAGE OF YEAR ATTENUATION EXCEEDS ORDINATE 










For the first case of distribution extension, Bayes' rule 
indicates how the extension is to be performed. Namely 


p ( R i > R 0 I T 1 T 2 “ P(T - 12 “ Tl) = 

P(T<T 2 -Tx i R, >R 0 ) > R 0 )> 


which states that the (probability that the instantaneous rain rate R i is 
greater than or equal to some specified rain rate R Q during the measurement 
time T less than or equal to the period from T x to T 2 ) times (the probability 
that T is less than or equal to T 2 - T J equals (the probability that a 
measurement was made given that R.. _> R Q ) times (the probability that R^ 2L R Q ). 
For the obvious limit that T = T 2 - T 1S this results in the identity 

p ( R i 1 R 0 I T = T > ' T 1 ) (1) = {1) P(R i - R o } ‘ 

The term P(T < T 2 - T a | R.. _> R Q ) accounts for the fact that if one made a 
nor.rep resen tati ve set of measurements during T when R.. _> R Q this term would 
correct the measured distribution to the true distribution. The desired 


quantity is 


p(R, > R 0 > 

V 


P( R i 1 R o ! T 1 T 2 - T i) P( T 1 T 2 - Tj) 
P(T£T 2 - T ! ! R-j > R 0 ) 


For the second situation, namely measurements were made during a 
time T between T x and T 2 and one wishes to extend the distribution to a longer 
sampling period T 0 to T 3 , where T 0 <^T X < T 2 .£T 3 , the same formula applies 
except Ti becomes T 0 and T 2 becomes T 3 . Here, for example, the term 
P(T < T 3 - T o | R i > R ) will account for the seasonal frequency dependence 
of thunderstorms, etc. The value of this term is difficult to determine, 
however the frequency of thunderstorms and number of hours with rain in excess 
of a -specified level will give approximations to the P(T < T 3 - T 0 | Rj > R Q )- 


Referring to Figure 3-24, the cumulative rain distribution has been 
determined (Step 18) by distribution extension. If the attenuation statistics 
were taken at the operating frequency, the cumulative attenuation statistics 
may be plotted. The construction of this curve, P(AL^ _> AL q | T < T 3 - T e ) , 
is shown in Figure 3-27. In this figure, the P(R i >_ R Q _ | T < T ? - T {) , 

P(R. > R 0 | T < T 3 - T 0 ) and P(AL i > AL Q | T < T 2 - T j)- are available from, the 


* Note: In this section l_ e , the effective path lengths, has been abbreviated L. 
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RAIN RATE, R Q 

i -- J X 1 1 1 I I — l 1 1 1 

RAIN ATTENUATION, AL Q 


FIGURE 3-27 , CONSTRUCTION OF CUMULATIVE ATTENUATION STATISTICS 
USING THE DISTRIBUTION EXTENSION TECHNIQUE 
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^measurements and the distribution extension. The desired curve 
P (AL^ _> AL q j T _< T 3 - T 0 ) is constructed by noting that the distance 
c - d is equal to a - b. For example* if the rain rate during the attenuation 
measurement period T exceeded R 1 mm/hr with probability P 1 (point a), and 
the measured attenuation^ for the same probability was about A' dB (point b) , 
and the rain rate exceeded R' for probability P" of the total period T 3 - r c 
(point c)> then the attenuation for the entire period T 0 to T 3 would be expected 
.to exceed A 1 dB for that same probability P" (point d). Using this technique, 
the entire curve P(AL 2 > AL 0 1 T 1 T a " T c ) generated. The probabilities P 
during the period T 0 to T 3 may also be expressed as a percentage of the time 
T 3 - T 0 or, e.g., as the number of minutes during the period T 3 - T 0 . 

3.5.4. 2 Frequency Scaling . If Steps 19B, C or D are the case, Steps 20, 21 
and/or 22 may be required. Step 20 is carried out by one of the following 
frequency scaling techniques. These are 

. 20-1) empirical scaling 

20-2) specific attenuation scaling 

20-3) Gaussian rain distribution scaling. 

Technique 20-1 utilizes experimental data to relate the total rain- 
, induced attenuation at one frequency to another frequency. The attenuation 
also depends on polarization because of the ellipsivity of the raindrops, but 
this effect is small compared to the accuracy of the above estimates and 
results primarily in depolarization. A general rule for scaling data based 
on 11 and 19 GHz radar data is (Ref. 3-24) 



where AL is the total attenuation along the path and f.. are the operating 
frequencies for the measurements. 
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Technique 20-2 utilizes the specific attenuation data (Ref. 3-4) 
in Section 3.1.3 to scale A(f 2 )/A(f 1 ). L is assumed to be independent of 
frequency. The result is 

• A(f?) _ a(f?) R (b(f 2 )-b(f 1 )) _ a(f 2 ) 

Wit attjy " attiT 

since b(f 2 ) - b(fi)~0. Note that this result is dependent on rain rate R. 


Finally, technique 20-3 assumes that each cell is spatially distributed 
along tl^e earth-space path according to the relation (Ref. 3-28) 


R(s) = R o exp(-s 2 /s.^) 


where s is along the path, z is the characteristic cell dimension and R 
the peak rain rate in the cell. Neglecting the term £ 
the result is 


i) 


is 





b(f? )-b(f i ) 
b(fi) 


which is imp! i city dependent on R through the term AL(f 1 ). 

Dropping the term 5> o "^2 ))/t>(Ti)^ can resu n j n errors when 
frequency scaling over wide frequency ranges. For a typical value of = 5 km, 
the error in scaling from 11.7 to 31.65 GHz is 20% due to this term alone. 
Therefore, for frequency scaling of more than two to one, it is recommended that 
.the full relation 


AL(f 2 ) 


a(f?) 



b(f 2 )/b(f i ) 

b(f 2 ) , 


Jh - b(f 2 )/2b(f 1 )) . 


t (1 - b(f 2 )/b(f 1 )) AL(fi) b(f 2 )/b(f 1 ) 

be utilized. 
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Tables of these frequency scaling factors are presented here for the 
common frequencies of 11.7, 15.3, 20, 28.56 and 31.65 GHz (see Tables 3-6 and 
3-7). The values of a{f 1 ), etc., are taken from the analytic relations in 
Section 3.1.3. 

3. 5. 4. 3 Elevatfon Angle Scaling . Step 21, the elevation angle scaling 
between the operational elevation angle 9 Q p and the measured data angle 9,- 
may be corrected for by using the relation 

AL(0 ) csc9„„ sinQ 

x op _ op _ 

AL (0) cscG sin9 

The corrections for the frequency and elevation angle scaling are first 
applied to the P(AL i > AL q | T < T 2 - T>) curve by varying the value of AL q , 
e.g., AL q becomes AL^ due to frequency scaling and elevation angle corrections. 
The curve P(AL.. > AL^ | T £ T 3 - T q } is then computed using the procedure 

^shown in Figure 3-24. 

! ? 

^ If more accurate cumulative attenuation statistics are required, it 

will be necessary to increase the attenuation measurement time T to obtain a 
truly representative sample of measurements. 

3. 5. 4. 4 Example of Attenuation Estimate Given Limited Rain Rate and 
Attenuati on Stati sti cs . The distribution extension technique has been applied 
to a six- month period from July to December 1974 at Rosman, NC. Here 
attenuation measurements were made for T - 608 minutes during the T = 725 
minutes of measurable rainfall. Rainfall measurements were then made from 
July (T ) to December (T 3 ). Using the technique in Figure 3-27, the 
attenuation distribution for the T 3 - T 0 period was estimated (see Figure 3-28). 

In this case steps 19 and 20 were not required, since the measurements 
were made at 20 GHz and at the ATS-6 'elevation angle. If an estimate were 
desired for another satellite at another frequency steps 19’ and 20 would be 
carried through using the attenuation statistics obtained from the distribution 
extension technique. 
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Frequency fj, GHz 


TABLI: 3-u 

FREQUENCY SCALING RATIOS FROM SPECIFIC ATTENUATION 
(Technique 20-2) 


Frequency f 2 > GHz 



11.7 

15.3 

20 

28.56 

31.65 

11.7 

1 

- 0 .024 

1.91R 

- 0.048 

3.66R 

- 0.007 

8.67R 

- 0 . ] 36 I 

11. OR 

15.3 

0.024 

0.52R 

1 

- 0.024 

1.91R 

-0 .083 

4.53R 

- 0 . 11 . 1 
5.76R ; 

20 

0 . 48 - 

0.27R 

0.024 

.52R 

1 

- 0.06 

2.37R 

- 0.065 ! 

3.01R : 

28.56 

.11 

0.12R 

0.083 

0.22R 

0.06 

0.42R 

1 

1 .27R _0 *°‘ 8 | 

31.65 

.136 

0.09R . 

.112 

0.17R 

0 .088 
.33R ' 

0 .028 

.79R 

1 1 


R in mm/hr 
f = 11.7 GHz 

a 

= 0.0162 dB*hr/km-mm 

b = 1.1641 

f - 15.3 

a 

= 0.0310 

b = 1.1401 

f = 20.0 

a 

= 0.0593 

b = 1.1165 

f = 28.56 

a 

= 0.1404 

b = 1.0568 

f = 31.65 

. a 

= 0.1787 

b = 1.0285 

a and b are derived 

from the analytic formulas in Table 3- 

3. More 

accurate values may 

be 

obtained by interpolating Table 3-2 

* 
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Frequency f l5 GHz 


TABLE 3-7 


FREQUENCY SCALING RATIOS FROM GAUSSIAN RAIN 


DISTRIBUTION SCALING 
{Technique 20-3) 

Frequency f z , GHz 



11.7 

15.3 

20.0 

28.56 

. 31.65 

11.7 

1.0 

1.858 

3.418 

7.408 

9.046 

15.3 

0.533 

1.0 

1.878 

4.198 

5.203 

20.0 

0.283 

0.533 

1.0 

2.328 

2.928 

28.56 

0.113 

0.218 

0.418 

1.0 

' 1.296 

31.65 

0.0833 

0.168 

0.318 

0.778 

1.0 


(b(f 2 )/b(f,)) 


Note: 3 = AL (f i ) 


and Zo = 4 km 





PERCENT OF TIME ABSCISSA IS EXCEEDED 


RA1NRATE (mm/hr) 

10 20 30 40 50 60 70 



FIGURE 3-23: ATTENUATION STATISTICS BASED ON DISTRIBUTION 

EXTENSION TECHNIQUE (FIGURE 3-24) . 
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4.0 CLOUD, FOG, SAND AND DUST ATTENUATION 


4.1 SPECIFIC ATTENUATION OF WATER DROPLETS 

The attenuation of clouds and fog arises from the small droplets (gen- 
erally less than 0.01 cm diameter) and can be computed in the Rayleigh approxi- 
mation, Therefore specific attenuation is related to the water vapor density 
of the cloud or fog (Ref. 4-1) as described in the Appendix. The resulting 
specific attenuation for clouds or fog is 

A c = K c % dB/km 

• where K is the specific attenuation per unit water vapor density p . K is 
c 3 3 w c 

in the units dB*m /g*km or dB*m /kg-km depending on the units of p w « 

. ' Values of K c versus frequency are given in Figure 4-1 (Ref. 4-2). .In 

engineering calculations for fog, which is generally most intense at the surface 
on a cool evening, the 10°C curve is recommended for summer and the 0°C is more 
appropriate for other seasons. For clouds, as with rain, the 0°C curve is- recoir- 
- mended. Note that these are theoretically-derived curves and so" their value is 
unverified. 

4.2 CLOUDS* 

O 

For heavy cumulus or cumulonimbus clouds, p is -typically between lgm/m 
3 w 3 

and 2.5 gm/m (Ref. 4-3),. although values as high as 4 to 5 gm/m are occasionally 

*Much of this work is taken from ORI Tech. Report 1187, D. Rogers and R. Kaul , 

June 1977. 
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SPECIFIC ATTENUATION K c (dB»m 3 fgm»km) 



FREQUENCY (GHz) 


IGURE 4-1. THEORETICAL ATTENUATION B? WATER CLOUD AT VARIOUS TEMPERATURES AS 
A -FUNCTION OF FREQUENCY 



observed (Ref. 4-4). Fair weather cumulus clouds have liquid water contents 

about an order of magnitude smaller (Ref. 4-5). For nonprecipitating events an 

3 

upper limit of 2.5 gm/m seems appropriate for system design in the U.S. and 
Canada. More detailed calculations may be made using the data in Table 4-1 
(adapted from Ref. 4-6) and local cloud data. 

Estimates of the path length L c through the clouds are more difficult 
to make. The maximum dimensions (horizontal and vertical) for heavy cumulus 
and cumulonimbus clouds are about 10 km (Ref. 4-7). If one assumes a right 
circular cylinder configuration and knowing the elevation angle, an approximate 
length l_ c may be determined. However for elevation angles below 30 degrees the 
effects of multiple clouds should be considered. As a general rule, a maximum 
L c * 10 km is not unreasonable, resulting in total cloud attenuation '(upper ^ 
bound) at 11.7 GHz of 


K c P v /tl max = (°- 13 d8m3 /km ■gm)(2.5 gm/m 3 ) (10 km) 

= 3.25 dB 

This result is quite high and should be consi dered an upper bound. The cor- 
responding upper bounds at 20 and 30 GHz are 10 dB and 20 dB, respectively. 

- Measurements made at Slough, United Kingdon, do not confirm these re- 
sults^(Ref. 4-1). The vertical attenuations at 95 and 150 GHz made between 
October 1975 and May 1976- are shown in Table 4-2. 


. TABLE 4-2 

ADDITIONAL ATTENUATION IN CLOUDS 


Jm. 

i Attenuation (dB) 

f=95 GHz 

f=150 GHz 

Strato-cumul us 

1 0.5 -1 

0.5 - I 

Small, fine weather cumulus 

0.5 

0.5 

Large cumulus 

1.5 

2.0 

Cumul o-nimbus 

2-7 

3-8 

Nimbus-strato (rain cloud) 

2-4 

' 5-7 ‘ 


Note: Water vapor density at ground level, was 4 to 11 gm/m 3 resulting in a 
clear-air attenuation of 0.7 to 1 dB and 1 to'3'dB at 95 and 150 GHz , 
respectively. 
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TABLE 4-1 

PROPERTIES OF CLOUDS 


Cloud Type 

Water 
Density 
( q/m 3 ) 

Base 

Height 

(m) 

Top 

Hei ght 
(R») 


Cirrostratus, Arctic 

0.1 

4000 

6000 


Cirrostratus, Midlatitude 

0.1 

5000 

7000 


Cirrostratus, -Tropical 

0.1 

6000 

8000 


Altocumuls and Altostratus 

0.15 

2400 

2900 


Low-lying Stratus 

0.25 

150 

1000 - 


Stratocumulus 

0.25 

1320 

660 


Fair Weather Cumulus 

' 0.5 

500 

1000 



1.0 

1000 

1500 



, 0.5 

1500 

2000 


Cumulus with Rain, 2.4 mm/hr 

0.1 

0 

500 



1.0 

500 

1000 


• 

2.0 

2000 

3000 


Cumulus with Rain, 12 mm/hr 

0.5 

0 

400 



2.0 

400 

1000 



4.0 

1000 

4000 


Cumulus Congest us 

0.3 

1000 

1200 



0.5 

1200 

1600 


• 

0.8 

1600 

2000 . 



1.0 

2000 

2500 



0.5 

2500 

3000 

— 

Cumulonimbus with Rain, 


i 



150 mm/ hr 

6.3 

0, 

300 



7.0 

300 

1000 


i 

8.0 

1000 

4000 


» 

4.0 

4000 . 

6000 



3.0 

6000 

8000 



0.2 

8000 

10000 - 
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These results, assuming the data in Figure 4-1 is correct and l_ c = 10 km, yield 
values for p = (8 dB ) / ( 5 dB- m 3 /kg-gm)(10 km) 

w o 

= 0.16 gm/m J . 

If one assumes p w = 2.5 gm/m 3 , then 

L c = (8 dB ) / ( 5 dB-m 3 /kg-gm) (2.5 gm/m 3 ) 

= 0.64 km. 

Neither result seems reasonable and so some engineering judgement must be used 
here for the' total cloud attenuation. 

Measurements made near Austin, Texas (Ref. 4-2) at frequencies of 35 
and 95. GHz indicated mean attenuations of 0.18 and 0.61 dB, respectively for 
strato- cumulus clouds. Correspondingly for cumulus clouds (probably ranging 
between the small and large cumulus values in Table 4-2), attenuations between 
0.12 and 0.34 dB were observed, respectively. The maximum cloud attenuations 
at 35 and 95 GHz were 2.3 and 10 dB, respectively, arising from rain-bearing 
clouds . 

Measurements of cloud attenuation can be most easily made with radio- 
meters and verified with satellite beacon measurements. These attenuation values 
fall in the range where radiometers, are most accurate. 

4.3 FOG 

Fog results from the consideration of atmospheric water vapor into.water 
droplets that remain suspended in air. When the resulting cloud of water droplets 
or ice crystals envelops an observer and restricts his ability to detect a black 
target against the sky horizon background by the 2 to 5 % threshold contrast of ■ 
the human eye at a distance of one kilometer or less, the international defini- 
tion of fog has been satisfied (Ref. 4-8 and 4-9). The relation (Ref. 4-10) 

p w = 0.0156 V -1 * 43 

appears to adequately relate p w to the visibility V in km for p w <_ 0.4 gm/m . 
This applies only to advective (coastal) fogs which form over open water as a 
result of the advection (movement) of warm, moist air over colder water. For 
radiation (inland) fog which forms in air that has been over land during daylight 




hours preceding tire night of formation, the relation (Ref. 4-11) 

P w - (42 V )' 1 - 54 

has been derived. The conditions favorable for radiative fog are: (1) air that 
has been under a cloud cover (with or without precipitation) during the day 
previous to its formation and (2) pools of air cooled to an excessive degree 
which collect in valleys, low, marshy land and along rivers on calm, clear 
nights. The maximum p w is generally less than 1.0 gm/m 3 . 

The attenuation due to fog is modelled by the same relation as clouds, but 
using the 0 or 10 C value for the specific attenuation per gm/m 0 . Unfortunately, 
the spatial distribution of p w or l_ c has not been measured via satellite beacons. 
Typical p w values for fog and haze are given in Table 4-3. 


TABLE 4-3 

PROPERTIES OF FOG AND HAZE 


Type 

Water 
Densi ty 
gm/m^ 

Base 
Hei ght 
(m) 

Top 

Height 

(m) 

Visibility 
(m) I 

Fog Layer* 

0.15 

0 

50 

i 

! 

Radiation (inland)Fog^ 

0.11 

- 

- 

100 - 

Advection (coastal) 

Fogt 

0.17 

. 


200 

Haze, Heavy* 

0.001 

0 

1500 

- 

*Adapted from Ref. 4-6. . 

tAdapted from Ref. 4-8. 
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4.4 


SAND AND DUST ATTENUATION 


Sand and dust scatter electromagnetic energy and their effect may be 
evaluated via Mie scattering. To date simulated measurements have been carried 
out in the laboratory (Ref. 4-12). At 10 GHz and concentrations of sand and 
dust less than 10~ 15 gm/cm 3 the measured specific attenuation was less than 
0.1 dB/km for sand and 0.4 dB/km for clay. Severe storms have concentrations . 
exceeding these values. 

Blowing sand and dust storms occur in some regions of the U.S. These 
are recorded by the Weather Service as part of the Local Climatological Data (LCD) 
at the 291 stations. Ground stations needing this information should review the 
data recorded by a nearby LCD recording station. 

y' 

The vertical extent of these sand storms is unknown, but it seems un- 
likely that high concentrations would exceed 1 km. The path length is expected 
to vary between h and 3 km, generally resulting in a total additional attenuation 
due tO’ sand of the order of 1 dB or less. No measured satellite beacon link data 
is available to confirm these results. 


Qp 
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APPENDIX 

CONVERSION OF WATER VAPOR PARAMETERS 

The U , S . and Canadian weather services provide measurements of the 
surface relative humidity. To calculate the attenuation due to water vapor, it 
is necessary to convert this quantity to the mass of water in a unit volume of 
air (gm/m 3 ) at the surface. In calculations of the refractivity the partial 
pressure (mbar) of water vapor is frequently required. In this appendix the 
conversion formulas are presented. 

CONVERSION OF RELATIVE HUMIDITY TO PARTIAL PRESSURE OF WATER VAPOR 

To a very good approximation the earth's atmosphere obeys an ideal gas 
law equation of state, except during periods when the water vapor will be super- 
saturated or supercooled. Fortunately, except for cloud physics, this is not a 
contributing factor for earth-space propagation. 

The equation of state for dry air is pv = R^T where R^ = 287 joule/kg • °K, 
p is the pressure in newtons/m 2 (1 newton/m 2 = 1.45 x 1 lO" 1 * lb/in 2 ), v is the 
specific volume in m 3 /kg ( v = 1/p, p = density) and T is the temperature in °K. 

i 

The equation of state for water vapor at atmospheric pressures and 
temperatures has been experimentally determined to be 

• ev * 461 T 

where e is the partial pressure of the water vapor. 


A-l 
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Since both dry air and water vapor act like perfect gasses, the 
Gibbs -Dalton law may be utilized for humid air so that 


H d R d + MvaA 
V M d + M w / 


T = (p d + e) v 


where p, v and T are the total pressure, specific volume and temperature of 
the humid air. and M are the masses of dry air and water vapor in grams, 
and = 287 and = 461 joule/kg * °K, respectively. 

The relative humidity is defined as 


where e g is the saturation vapor pressure corresponding to the temperature of 
the mixture. The saturated partial pressure may be found in the ASME Steam 
Tables or the Handbook of Chemistry and Physics. A plot of e $ in newtons/m^ 
versus temperature is given in Figure A-l. 

CONVERSION OF RELATIVE HUMIDITY TO WATER VAPOR DENSITY 

The water vapor density P w (kg/m 3 ) at a given temperature T and vapor 
pressure e may be found from the ideal gas law; i.e., 

» w ■ ■ «/(V> = < R - H -)e s /(y) 

where e is the vapor pressure in newtons/m 2 , R w = 461 joule/kg • °K, and T is the 
absolute temperature. 

Alternatively, the saturated water vapor specific volume v w is given 
in the ASME Steam Tables and the Handbook of Chemistry and Physics. The con- 
version for the specific volume is 

vjm 3 /kg) = 0.062 v (ft 3 /lbm) 
w V/ 

where Ibm is the mass equivalent to a pound in a standard accele- 

ration field. Then p w = (v J" 1 . 
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FIGURE A-l . THE SATURATED PARTIAL PRESSURE OF WATER VAPOR VERSUS TEMPERATURE 
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^APPROXIMATE CONVERSION OF RELATIVE HUMIDITY TO COLUMNAR WATER CONTENT 

For radiometric measurements, it is useful to know the total amount 
of water vapor in a column along the line of sight. To a first approximation, 
one can assume that the water vapor density p w is distributed exponentially 
from the surface to the 0°C altitude. The zonally averaged height of this 
isotherm is given in Figure A-2 for the rain climates shown in Chapter II. 

The amount of water in a column along the line of sight is 

A 

R = cscO p w f e ~ t(x 'A 1 ) dx 

5 A X 


where Ai N is .the ground station altitude, A,, is the 0°C altitude, p w - (R.M.) e $s 
is the surface water vapor density and y is the scale height. The y is deter- 


mined by the boundary conditions 


Y = - 


Inp, , - Inp,, n 
w s w Q o c 

\ " A 2 




where p ' and p are computed from Figure A-l for the appropriate 


' w s """ W 0°C 


surface temperature. The total column of water is then 

{cscG)(p - Pw 0 o c H a 2 - A X) 

R kg/m 2 = -y- — T } ■ 

In p - In p 
w w Q o c 


AN EXAMPLE 

The following example applies to the Rosman, N.C. ground station 
(elevation = 880m, latitude = 35°N) to ATS-6 (elevation angle = 47°) for the' 
month of September (T surface = 80°F = 26.7 °C, R.H. = 60%). 

i 

At this temperature the saturated water vapor content e g is 0.51 psia 
or 3516 newtons/m 2 (from Figure'A-1), so the surface vapor pressure is 2110 
newtons/m 2 . 
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MODEL FOR RAIN RATE CLIMATES A THROUGH H 
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SEASONAL 1 - WINTER (NORTHERN HEMISPHERE) 

2 - SPRING 

3 - SUMMER 

4 - FALL 

FIGURE A-2. LATITUDE DEPEUDENCE OF ZO'iALLY AVERAGED 0°C ISOTHERM l.EI 
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The water vapor densities are computed from the. ideal gas law at 
the ground station altitude Aj . 

(R.H.) e. 


pw s * R 


w 


•• fcj } w * 1 - 53x10 ' 1 k9/m3 


' w °° c - = 2 - 9 * io " 3 k9/m3 

The altitude of the 0 °C isotherm is approximated from Figure A-2 
A 2 - 3700 m 

Using the formula to calculate the columnar water content yields 
R=csc(47°) ^n°i 5 8 3 ^in '2 9 2 ' 9 ^ X 10 - 3 = 28 * 75 kg/m2 
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1.0 PATH DIVERSITY 


1.1 THE PATH DIVERSITY CONCEPT 

Rain attenuation often degrades earth-space paths operating above 
10 GHz so seriously that the requirements of economical design and reliable 
performance cannot be achieved simultaneously. To overcome this problem, Hogg 
(1968) proposed the use of path diversity on earth-space paths to achieve the 
desired level of system reliability at a reasonable cost compromise. This 
proposal was based on the hypothesis that rain cells and, in particular, the 
intense rain cells that cause the most severe fading are rather limited in 
spatial extent. Furthermore, these rain cells do not occur immediately 
adjacent to one another. Thus, one might expect that the probability of simul- 
taneous fading on two paths to spatially separated earth terminals would be 
less than that associated with either individual path. This hypothesis was 
tested first by Wilson (1970) using-radiometric noise emission measurements to 
determine the rain attenuation on separated paths and by Hodge (1974a) using 
actual earth-space paths. These and other ensuing experiments have clearly 
demonstrated that path diversity is an effective technique for improving system 
reliability in the presence of rain attenuation. 

The reliability improvement resulting from the use of path diver- 
sity must ultimately be balanced against the increased system costs arising 
from the need for a second earth terminal site, antenna, equipment, and 
connecting link between the two earth terminals. Of course, redundant equip- 
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merit normally utilized at a single terminal might simply be separated physically 
to achieve this improvement at a reduced cost. Also* a judicious design might 
incorporate the connecting link as part of the trunk link to the earth terminals 
in order to further reduce the overall cost. Finally, the overall system 
reliability degradation due to the requirement of additional system equipment 
must be weighed also. This cost and equipment reliability will not be con- 
sidered further in this discussion. The remainder of this discussion will bq 
directed toward the design and performance of path diversity systems independent 
of cost and equipment reliability aspects. 

A typical path diversity earth terminal configuration is shown in 
Figure 1.1-1 along with the definitions which will be used in the following dis- 
cussion. The parameters are defined as follows: 

Az = azimuth of earth-space path (degrees) 

EL = elevation of earth-space path (degrees) 

d = separation between earth terminals (km) 

3 = orientation of earth terminal baseline (degrees 0<& < 180) 

-0 = major axis of predominant rain cell orientation degrees, 

0 _< 0 < 180) • 

Two alternative approaches (Brandi nger-1978) have also been suggested 
•as techniques for improving reliability in the presence of rain attenuation; 
these are angle and frequency diversity. The first, angle diversity, uses one 
ground site with more than one earth-space path to satellites located in sepa- 
rated orbital positions; thus, the paths are oriented along different azimuths. 

If a 'rain cell is located on one path at some distance from the terminal, the 
result will be quite similar to that for path diversity; however, if a rain cell 
is located near the terminal, little improvement results. This approach is not as 
effective as path diversity but may find utility in cases where multiple satellites 
are available. 

t 

The second approach, frequency diversity, is based on the frequency 
dependence of rainfall attenuation. In this case channel 'assignments both above 
and below 10 GHz are incorporated in the same system. Thus, in the event of rain- 
fall, high priority traffic is diverted to the lower frequency which is less sus- 
ceptible to rain attenuation. Thus, the channel capacity is adaptively reduced to 
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Figure 1.1-1 Path Diversity Configuration 
and Parameters 
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maintain reliability for a portion of the traffic. This is equivalent to 
assigning all high priority traffic to a lower frequency and tolerating a 
lower reliability channel above 10 GHz. This approach will probably not be 
useful in general applications. 

1.2 DIVERSITY GAIN AND DIVERSITY ADVANTAGE 

. In order to characterize the performance of a path diversity system, 
it is convenient to establish a descriptive parameter for -this purpose. Two 
such parameters have been proposed and utilized in the literature; they are 
diversity gain, G, and diversity improvement, I. 

Let us consider the distribution of rain attenuation, i.e., excess 

path attenuation, fade depth, or attenuation above the clear sky level, exceeded 

*f* 

on a single path, A(T), for a given percentage of time, T, and the rain 
attenuation, A^. (T) ,' exceeded jointly on two separated paths for the- .same per- 
centage of time as shown in Figure 1.2-1. It appears to be more customary in 
the literature to plot- A along the abscissa and T along the ordinate, even 
though A is generally taken to be a function dependent upon T. Obviously, if 
the attenuation distribution is measured or modelled properly, the functional 
.relationship between A and T in any given experiment is unique whether one 
chooses to think of A being dependent upon T, A(T), or T being dependent upon 
A, T(A) . Furthermore, given A(T), T( A) is immediately known. Therefore, both 
functional relations will be utilized in the following. 

■ It should also be noted that percentage time , T{A), will be assumed 
to be equivalent to the probability- of occurrence, P(A), multiplied by 100, T(A) 
= 100 P( A) . This, of course, is only an approximation for experimental data of 
limited duration. Therefore, the term percentage time will be utilized in the 
following discussion instead of probability in order to emphasize this fact. 

Diversity gain, G, (Hodge-1974a) may be defined as the difference 
between the rain attenuation exceeded on a single path and that exceeded 
jointly on separated paths for a given percentage of time, i.e., 

G(A) = A(T) - A div (T) (1.2-1) 


In this section A denotes the total attenuation. AL was used in 
section II. • e 
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^Diversity advantage, I, (Wilson and Mammel-1973) may be defined as the ratio of 
the percentage of time exceeded on a single path to that exceeded jointly on 
separated paths for a given rain attenuation level, i.e., 

T( A) 

1(A) = (1.2-2) 

W A > 

Diversity gain may be interpreted as the reduction in the required system margin 
at a particular percentage of time afforded by the use of path diversity. Al- 
ternatively, diversity advantage may be interpreted as the factor by which the 
fade time is improved at a particular attenuation level due to the use of path 
diversity. 

Clearly, both diversity gain and advantage are functions of attenua- 
tion level or, implicitly, percentage of time. Furthermore, referring to Fig. 
1.2-1, it is evident that the specification of either diversity gain or improve- 
ment along with single terminal fade distribution permits the joint attenuation 
|>to be reconstructed. Therefore, in principle, these two parameters are equiva- 
lent and contain the same information; however, in actual practice this is not 
the case. This can be seen by again referring to Fig. 1.2-1 and also observing 
that the uncertainty of experimental data increases as the percentage time of 
occurrence decreases. , Now, the diversity gain determined at a particular fade 
depth depends upon two quantities, A(T) and A di - V (T), subject to uncertainties 
of the same order. In contrast, the diversity advantage associated with the 
same fade depth. A, will be a ratio of a quantity, T(A), subject to the same- 
order of uncertainty as that determining G d ^ v (A), to a quantity, T div (A), 
subject to a much larger uncertainty. Therefore, one may conclude that the 
diversity advantage factor will be subject to greater uncertainty than the 
diversity gain parameter associated with the same fade depth. The variability 
of the diversity advantage factors as. compared with the diversity gains reported 
in the experimental literature and referred to in the following readily confirms 
this argument. 
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In the preceding discussion, it has been implicitly assumed that 
the fade distributions associated with each individual earth terminal were 
identical. However, in actual practice this is seldom the case. These distri- 
butions may differ because the measurement period is not long enough for the 
spatial rainfall distribution to become uniform; in this case the differences 
between the distributions are an indication of the uncertainty of the measure- 
ment. Alternatively, local climatological variations may give rise to real 
differences between distributions measured only a few kilometers apart; very 
little is known definitively about this latter effect. In general, however, 
since there is little or no basis to give greater weight to one distribution 
over another, the definitions given in eqs. 1.2-1 and 1.2-2 should be geneia;- 
ized for experimental data analysis to: 

G(A) ■ A ave (T) - A aiy <T) U-3i 


1(A) ■ 


^ave 

T div 


(A) 

W 


O /i\ 


where A (T) is the average of the single terminal fade depth exceeded for 
the percentage time, T, and T ave (A) is the average of single terminal percent- 
ages of time exceeded at the fade level A. 

.1.3 DIVERSITY EXPERIMENTS 

1.3.1 Experimental Methods 

Four methods have been used to experimentally determine the per- . 
formance of earth-space path diversity systems. The first, of course, is the 
direct measurement of path attenuation on actual earth-space paths using 
synchronous satellite beacons. This method is the most direct method but is 
expensive and relies on the availability of appropriate satellite beacons. 

The second method simulates path attenuation from measurements 'of radiometric 
noise emission. This radiometric technique is based on- the simple relationship 
between radiometric sky temperature and path attenuation which is applicable 
when scattering is not significant, i.e., in rainfall at frequencies below 
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about 30 GHz. This method is generally less expensive than direct satellite 
Deacons. However, as the path attenuation increases, the sky temperature 
asymptotically approaches the mean absorption temperature of the path, about 
270°K; this limits the useful dynamic range of the inferred attenuation levels 
to approximately 10-12 db. 

The third technique uses conventional meteorological radar measure- 
ments of backscatter from rainfall to estimate in turn the reflectivity, rain 
rate, and specific attenuation associated with each resolution cell. Path 
attenuations may then be found by summing the specific attenuation along the 
paths of interest. This approach is subject to uncertainties resulting from 
variations in the drop size distribution and the presence of ice, hail, and 
snow on the path. Nevertheless, the number of possible diversity configurations 

s' 

and paths that can be synthesized from digitally recorded reflectivity data makes 
this technique extremely attractive. It appears to be the only practical 
method for the study of spatial anisotropy, local climatological, and baseline 
orientation questions. 

The fourth method utilizes radiometric measurements of solar emission 
} to determine the atmospheric path attenuation. This method does not suffer from 
the dynamic range limitation of the radiometric method mentioned above. However, 
by virtue of the suntracking, the measurement paths are limited to a unique rela- 
tionship between the path azimuth, path elevation, and local time of day. This 
makes the separation of dependences on these variables virtually impossible. 
Therefore, this latter technique will not be considered in detail in the following. 

1 . 3 .-2 Experimental Results 

The experimental diversity results available in the literature are 
summarized in Table 1.3. 2-1. This table includes the results reported for each 
of the four methods described in the preceding section. In each case the 
reference is cited along with the location of the experiment, the frequency, 
separation distance, baseline orientation, path azimuth, and path elevation. In 
cases where multiple measurements are reported, the range- of the appropriate 
parameters is indicated. 
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' TABLE 1.3. 2-1 SUMMARY OF DIVERSITY EXPERIMENTS 


I. SATELLITE EXPERIMENTS 


REFERENCE 

LOCATION 

FREQ. 

SEPARATION 

BASELINE 

Az 

Ik 

Hodge (1974) 

Columbus, Ohio 

15.3 GHz 

4. 0-8.3 km 

159-164° 

210° 

38° 

Hodge (1976b) 

Columbus, Ohio 

20-30 

13.2-14,0 

33-151° 

197° 

40° 

Vogel, et al (1976) 

Austin, Texas 

30 

11.0 

0° 

.172° 

55° 

Hyde (1976) 

Boston, Mass. 

18 

6.7-35.2 

74-93° 

212° 

36° 


Columbus, Ohio 

18 

5.1-38.9 

91-95° 

196° 

42° ■ 


Starkville, Miss. 

18 

8.3-40.0 

105°- 113° 

190° 

,„51° 


II. RADIOMETER EXPERIMENTS 




Wilson (1970) 

Crawford Hill, N.O. 

16 

3.2-14,4 

135° 

226° 

32° 

Wilson & Mammel (1973) 

Crawford Hill, N.O. 

16 

11.2-30.4 

135° 

226° 

32° 

Gray (1973) 

Crawford Hill, N.J. 

16 

19.0-33.0 

45-135° 

226° 

32° 

Funakawa & Otsu (1974) 

Kokubunji , Oapan 

35 

15.0 

— 

180° 

45° 

Hall & Allnutt (1975) 

Slough, England 

11.6 

1.7-23.6° 

20-106° 

198° 

30° 

Allnutt (1976) 

Slough, England 

11.6 

1.7-23.6° 

20-106° 

198° 

30° 

Strickland (1977) 

Quebec, Canada 

13 

18.0 

11° 

122° 

19° 


Ontario, Canada 

13 

21.6 

1° 

116° 

16° 

Bergmajin (1977) ' 

Atlanta, Georgia 

17.8 

15.8-46.9 

141-146° 

228° 

38° 


Denver, Colorado 

17.8 

33.1 

86° . 

197° 

43° 


III. 

RADAR EXPERIMENTS 




Goldhlrsh & Robison (1975) 

Wallops Island, Va. 

13-18 

2-20 

0-180° 

- 0-360° 

45° 

Goldhirsh (1975) 

Wallops Island, Va. 

13-100 

2-20 

0-180° 

0-360° 

45° 

Goldhlrsh (1976) 

Wallops Island, Va. 

18 

.2^20 

0-180° 

0-360° 

45° 

Hodge (1978) 

Montreal , Quebec 

13 

4-42 

0-180° 

122-240° 

19-40° 


IV. SUNTRACKER EXPERIMENTS 




Wolfsburg (1973) 

Boston, Mass. 

35 

11.2 

158° 

--- 

— 

Funakawa l Otsu (1974) 

Kokubunji , Oapan 

35 

15.0 

— 

— 

— 

Davies S Croom (1974) 

Slough, England 

37 

10.3 

67° 

--- 

— 

Davies (1976) 

Slough, England 

37 

10.3-18.0 

67-110° 

— 

— 
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1.4 


PATH DIVERSITY DESIGN FACTORS 


1.4.1 Separation Distance 

Diversity gain depends strongly upon the earth terminal separation 
distance, d. The diversity gain increases rapidly as d is increased over a small 
separation distance, i.e., up to about 10 km; thereafter the gain increases more 
slowly until a maximum value is reached, usually between about 10 and 30 km 
(Fig. 1.4. 1-1). This maximum value is generally quite close to that value 
associated with uncorrelated fading at the individual earth terminals. The 
maximum value has been referred to in the literature as an optimum value in the 
sense that it is the largest observed value. It may be argued that negatively 
correlated rainfall along the separated paths could produce higher values of 
gain; however, no cases of this occurrence have been reported. 


The uncorrelated diversity gain may be calculated from the single 
terminal attenuation distribution in the following manner. First, the condi- 
tional single terminal attenuation distribution, P c (A) s is formed given that 
fading is occurring on the path. This is necessitated by virtue of the fact 
that the absence of fading on separated paths is well correlated during clear 
weather (Goldhirsh-1976). The joint conditional probability of uncorrelated 

attenuation on separated paths, P , . (A), is then simply the square of the 

1v c 

single path conditional probability, i.e., 




(1. 4-1-1) 


or 



T C (A) 


7ioo 


(1.4. 1-2) 


where the relationship between probability and percentage time is as described 
above.- The relationships between the absolute and conditional quantities are 


P(A) = P r P c (A) 


(1.4. 1-3) 


or 


T (A) = t r t c (a)/ioo 


(1. 4.1-4) 
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DIVERSITY GAIN, G (dB) 



SITE SEPARATION DISTANCE, <1 (km) 


Eigure 1.4. 1-1 Diversity Gain, G, versus Separation Distance, d, 
f = 18 GHz. (Horizontal Dashed Lines Represent 
Opti.TVjm Labels) (Coi dhi rsii and Robison - I j 75 j 
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where the subscript R is associated with the probability or percentage time 
that rainfall attenuation is occurring on a single path. 

A practical difficulty is associated with the determination of 
P R and T r . Namely, as a result of clear air and/or cloud induced scintillations 
on the received signal it is usually not possible to clearly distinguish 
between small rain induced fades and clear air scintillation. Therefore, it 
has been customary to establish a threshold of attenuation beyond which the 
fading is assumed to be due to rain. This level has been typically taken to be- 
tween 0.5 and 3.0 dB. 

In contrast to the uncor.related case, one may argue that correlated 
fading may occur for paths separated by distances associated with typical rain 
cell separation distances. This would tend to decrease the diversity gain over 
some range of separation distances. Such an effect may be inferred from the 
rainfall statistics of Freeny and Gabbe (1969); however, these statistics are- 
associated with point rainfall rates rather’ than path average rainfall rates. 
Again, no definitive report of this effect has been published to date. 

1.4.2 Baseline Orientation 

The perpendicular separation between parallel paths is greatest 
when the earth terminals are located on a baseline perpendicular to the' pro- 
jections of the paths on the earth's surface. This arrangement minimizes the 
possibility of both paths passing through the same rain cell. Nevertheless, 
the dependence of diversity gain on baseline orientation is quite weak except, 
possibly, for very short separation distances (Fig. 1.4. 2-1). 

The baseline orientation problem is further complicated if spatial 
anisotropy of the rain cells, i.e., a preferred direction of rain cell 
elongation, is known to exist in the region of interest. In this case, a 
baseline orientation perpendicular to the preferred axis of rain cell orientation 
would be desirable if the direction of the propagation path were ignored. 

Considering both factors together, it appears that the most desirable 
baseline orientation is that which bisects the larger of the two angles between 
the projection of the propagation path and the preferred axis of rain cell 
orientation. 
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AVERAGE DIVERSITY GAIN (%) 


CASE I 


CASE !! 


CELL A CELL 





BASELINE ORIENTATION 
(DEG) 


BASELINE ORIENTATION 
(DEG) 


CASE ill 


CASE IV 




(DEG) 


CELL 




BASELINE ORIENTATION 
(DEG) 


Figure 1.4. 2-1 
Baseline 


Average Diversity Gain versus 
Orientation (Hodge-1978) 
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1.4.3 


Path Elevation Angle 


The separation distance required to achieve a given level of 
diversity gain increases as the path elevation angle decreases (Hodge-1978). 

This is due to the increased likelihood of path intersections with rain cells 
at lower elevation angles. This effect is coupled to the problem of rain cell 
anisotropy and path azimuth as noted below. 

1.4.4 • Path Azimuth Angle 

For synchronous satellites the path azimuth and elevation angles 
are not independent, and, thus, the dependence of diversity performance on these 
variables cannot be fully separated. If all rain cells were isotropic, one would 
expect no variation in diversity performance with azimuth angle other than that 
associated with the elevation angles. However, when rain cell anisotropy is con- 
sidered, there appears to be a weak improvement in diversity performance for 
path azimuths in the quadrant not containing the preferred axis of rain cell 

orientation. 

1.4.5 Link Frequency 

Experimental measurements to data have shown no significant depend- 
ency of diversity gain on the link frequency (Goldhirsh and Robinson-1975) over 
the 10-30 GHz frequency. Certainly the probability of a given attenuation level 
being exceeded on a single path is strongly frequency dependent. However, diver- 
sity gain is a conditional statistic based upon the difference between single 
terminal and diversity attenuation levels. However, for link frequencies above 
30 GHz, attenuation on both paths simultaneously can be sufficient to create ' 
an outage. Therefore, extrapolation beyond 30 GHz is not recommended. _ 

1.4.6 Anisotropy of Rain Cells Along a Front 

j 

There is a tendency for convective rain cells associated with 
frontal activity to occur in bands nearly perpendicular to the direction of 
movement of the front. The direction of motion of the cells within such a 
band tends to be along or slightly ahead of the direction of the front, and, 
furthermore, the more intense cells tend to elongate in~ their direction of 
motion (Harrold and Austin-1974). Thus, two types of anisotropy are evident. 
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The first is associated with the elongation of individual cells and is related 
to the probability of parallel paths passing through the same cells.- The second 
is associated with the statistics of the vector separation between rain cells 
and is associated with the probability of parallel paths simultaneously inter- 
secting two different rain cells. 

Fortunately, for simplicity's sake, these two preferred orientations 
are nearly parallel, and thus the same corrective action is required in each 
case. Namely, the baseline orientation should be nearly perpendicular to these 
preferred directions. 

It should be noted at this point that diversity gain is not 
directly dependent upon either the speed or direction of rain cell motion; it 
is dependent only on motion through the relationship between motion and spatial 
anisotropy. On the other hand, fade duration statistics will be quite dependent 
upon the speed and direction of rain cell motion. 

1.4.7 Local Climatology 

To a first order of approximation it is commonly assumed that the 
probabilities of rain cell occurrence are uniformly distributed over rather large 
regions of the earth's surface. This assumption may be invalidated by the 
presence of any one of the following features: mountains, large valleys, 

large bodies of water, or urban heat "islands". These features can give rise 
to nonuniform spatial distributions of rain cell probabilities. 

Spatial distributions of rainfall accumulation are readily available 
in the meteorological literature-; however, it is not currently known whether 
the use of these data is applicable to the question of earth terminal siting. 

For example, it may be argued that these rainfall accumulations are dominated 
by low rainfall rates and thus do not reflect the spatial distributions of 
intense rain cells which dominate the occurrence of high attenuation levels on 
earth-space paths. 

1.4.8 Switching Rates 

The fading rate on a single path is relatively slow. The highest 
rates reported are on the order of 0.1 dB/sec. at 15 GHz (Hodge-1 974a, Strick- 
land-1977). This implies that the decision and switching process for diversity 
paths may be quite slow,and should pose no significant problem in the system design. 

F 
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1.4.9 


Connecting Link 


The implementation of a path diversity system must incorporate 
a connecting link between the two earth terminals. If this link is closed, 
i.e., waveguide, coax, etc., its performances will be independent of meteoro- 
logical variables and will not directly influence the reliability improvement 
provided by' the use of path diversity. If, however, the connecting link 
'operates above 10 GHz in the atmosphere, the joint fading statistics of the 
connecting link with the earth-space paths must be considered. This degrading 
effect appears to be small except for cases of very long baselines or baseline 
orientations parallel to the earth-space propagation paths. (Ferguson and 
Rogers-1978) 

1.4.10 Multiple Earth Terminals 

Substantial link reliability improvements result from the use of two 
earth-space propagation paths. Thus one may conjecture that further improvement 
might result from the addition of additional diversity paths. Cost considerations 
generally make this approach unattractive. Furthermore, the resulting additional 
improvement is generally quite small. This can be' observed by considering the 
uncorrelated fading case discussed earlier. The conditional probability of 
joint fading on N paths, given that fading occurs on a single path, is 

P div = P >> (1.4.10-1) 

C N 

Determination of diversity gain for N diversity terminals shows that most of the 
gain is realized for two terminals with very little further increase in gain for 
additional terminals (Hodge-1974b) . 1 
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2.0 PATH DIVERSITY DESIGN INFORMATION 


2.1 EMPIRICAL MODEL 

The data available from early diversity experiments in New Jersey 
and Ohio (Hodge- 1974a, Wilson-1970, Wilson and Mammel-1973, Gray- 1973) were 
used to develop an empirical model for the dependence of diversity gain on „ 
separation distance, d, and single site attenuation, A. (Hodge-1976a) . The 
resulting model was of the form 

G=a(l-e" bd ) (2.1-1) 

where the coefficients a and b depended upon the single site attenuation 
according to 

a = A -3.6 (l-e"°- 24A ) (2.1-2) 

b = 0.46 (l-e"°- 26A ) (2.1-3) 

These results were nominally applicable to baseline orientations of 135°, i.e., 
NW-SE, elevation angles of 35°, azimuth angles of 220°, and link frequencies of 
15 GHz. The model reproduced all of the available data at the time within 0.7 
dB as shown in Fig. 2'. 1-1. The variation of the coefficients a and b are shown 
in Fig. 2.1-2. 

The coefficient a represents the diversity gain achieved for large 
separation distances according to Eq. 2.1-1. From Eq. 2.1-2, it is seen that 
this achieved diversity gain is simply A-3.6 dB for large fade depths. In 
addition, referring to Eq. 2.1-1, the quantity 3/b is. the separation distance 
required to achieve 95% of the maximum diversity gain for large fade depths; 
from the equation for b this model yielded a value of 6.5 km. 

! *The empirical relation should not be utilized for attenuation in 
excess of 15 or 20 dB, since it appears to reach an unphysical limit for large 
attenuation values such as those associated with frequencies above 30 GHz. 


*This and the following paragraph were written by R : Kaul , ORI , Inc. 
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d, SEPARATION (km) 


Figure 2.1-1 Diversity Gain versus Separation Distance 






Consider the following relation for the link margin: 


M' - M - A + G 


(2.1-4) 


where M' is the link margin at any time and M is the zero-rain link margin. 
In the limit where A becomes large 

M'*M - 3.6 - (A-3.6)e"°* 46d ^ -,_ 5 j 

-0.24A 

neglecting terms similar to e , etc. Thus for even moderate separation 

distances, d j> 10 km 

M' ~ M - 3.6 (2.1-6) 

which says that the link only degrades 3.6 dB for any level af attenuation. 

Rewriting Eq. 2.1-5 for A yields the amount of attenuation on 
the reference link allowable for a given M and M‘ , 

A = 3.6 + (M-M'-3.6)e +0,46d (2.1-7) 


Solving this for M = 20 dB, M' - n = 10dB for a low error rate and d = 10 km, 
yields 

A = 3.6 + (20-1 0-3. 6)e 0,46d (2.1-8) 

max 

~ 640 dB! 

Clearly this is a physically unrealizable number. A term must be added to Eq. 
2.1-7 and then into 2.1-1 to avoid this situation. The form of this term is 
now under study. 


A more recent analysis of radar data near Montreal, Quebec, 
(Hodge-1978) has yielded considerably larger values of 3/b; they range from_ 
13.8 to 31.3 km depending on the path azimuth and elevation angles. It is 
believed that these larger values are a result of the averaging method used 
to obtain the diversity gain in this case. The average percentage diversity 
gain was defined in this work to be 


ave 


t N G(A.) 
- A v 
N 


i=l 


A i 


x 100 


(2.1-9) 


and the coefficient resulting from the modeling of G avg = a(l-e ^ d ) are given 
in Table 2.1-1 . 
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TABLE 2.1-1 COEFFICIENTS FOR EMPIRICAL MODEL BASED ON MONTREAL DATA 

(f = 13 GHz, Hodge- 1978) 


CASE 

A z 

El 

a 

lb 

3/b 

I 

121.5° 

18.5° 

53.0% 

0.096 km -1 

31.3 km 

II 

181.5 

37.7 

51.2 

0.155 

19.4 

III 

240.1 

17.9 

50.5 

0.126 

23.8 

IV 

210.0 

40.0 

51.6 

0.156 

19.2 

2.2 

USE OF THE 

EMPIRICAL MODEL 





The experimental results in Figure 2,1-1 do demonstrate that the 
diversity gain does appear to apply up to the 8 to 10 dB level in the 10 to 
30 GHz frequency range. Using this fact, the cumulative attenuation statistics 
for a hypothetical ground station system of two identical stations separated 
by 10 to 15 km at Rosman, N.C., would yield statistics significantly better 
than one station. These are shown in Figure 2.2-1 based on the attenuation 
data obtained by the distribution extension technique in Section II. The 
diversity gain is the difference in attenuation levels for a given percentage 
of exceedance. For example, at A = 10 dB the diversity gain is approximately 
6.7 dB (see Figure 2.1-1). The diversity advantage (defined in Figure 1 .2-1) 
-would be about 4.8 at A = 3.3 dB. 
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FIGURE 2.2-1. PATH DIVERSITY GAIN STATISTICS FOr'rOSMAN, N.C. 
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1.0 INTRODUCTION 


1.1 MECHANISMS FOR SIGNAL FLUCTUATIONS AND LOW ANGLE FADING 

The amplitude, phase, and angle-of-arrival of a microwave signal 
passing through the troposphere vary due to inhomogeneities in the refractivity 
{clear air), clouds and ra.in. The effects occur, on time scales shorter than 
a minute and on spatial scales shorter than a kilometer. 

1.2 ANTENNA APERTURE EFFECTS 

The effects of the tropospheric turbulences and the antenna can not 
be totally decoupled because, of course, the measurements and operating sys- 
tems utilize antennas. The antenna aperture processes the incident wavefront 
with its spatial and temporal fluctuations into a received signal with only 
temporal variations. 

Wavefront tilt due to inhomogeneities and gradients in the refrac- 
tivity appear to the antenna as -an angle-of-arrival variation. Average- ele- 
vation angle ray bending is usually 10 times more pronounced than azimuthal 
ray bending. However wave tilt fluctuations tend to be randomly distributed 
in angle relative to the slant path porpagation direction, at least when 
the majority of the path is above the regime of surface- effects (surface ef- 
fects extend upwards several hundred meters). 
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Fluctuations occurring on spatial scales smaller than the size of 
the aperture are often referred to as wavefront ripple. This phase incoher- 
ence results in an instantaneous gain loss or degradation. 

The fluctuations described herein apply to the ground station downlink 
because its antenna is in close proximity to the turbulent medium. An uplink 
satellite path will suffer fluctuation gain degradation only due to scattering 
of energy out of the path. Because of the large distance traversed by the wave 
since leaving the troposphere, the wave arrives at the satellite antenna as 
a plane wave (no ripple) and with only minute angle-of-arri val effects. 
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2.0 AMPLITUDE FLUCTUATIONS BASED UPON THE -COMBINED 
EFFECTS OF AMPLITUDE AND ANGLE-OF-ARRIVAL 
FLUCTUATIONS — MODEL AND EXPERIMENTAL RESULTS 


2.1 OVERVIEW 

The phenomena of amplitude and angl e-of-arrival fluctuations combine 
to form received signal amplitude fluctuations. For many cases of propagation 
one or more of these effects may often be neglected. For example, a receiving 
system which employs an antenna with a wide beamwidth will not experience 
angl e-of-arrival -induced amplitude fluctuations for most elevation angles. 
Howeve/, such simplification is not always possible. The theory of wave propa- 
gation and scattering in random media allows a combination of the turbulence- 
induced effects to be performed in the context of weak fluctuation along a 
line-of-sight path. The work of Isbimaru (1978), which defines coherent and' 
incoherent field components as a plane wave propagates through a random medium, 
provides a method of combining amplitude and angl e-of-arrival effects into a 
model of received signal amplitude fluctuation. A model utilizing the concept 
of incident plane wave decomposition (see Figure 1) has been proposed by Theobold 
and Hodge (1978). 

2.2 VARIANCE OF RECEIVED SIGNAL AMPLITUDE 

The assumption of weak turbulence, i.e., log-amplitude variance 
ox 2 < 0.5, is invoked for a plane wave incident on a region of turbulence, 
propagating a distance L(km) and impinging on a circular aperture of diameter 
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d (meters). The antenna is assumed to have a Gaussian pattern function with 
half-power beamwidth B (degrees). If v is the received signal amplitude, an 
expression for signal variance relative to average power is 


S 2 = 10 loqi 


(v 2 / - ^v) 

<v > 2 


l c Jl + 8a|"iF(2) +~B r “ X i (4 a| ln(2) + B 


r c. + ! i 14 af In (2) + B‘ 


where 

I n - = 1 - exp [-L/Lo 
I Q = (1 1^) / (1 + af) 

ai 2 = amplitude variance 2 

A 

aa 2 = angle-of-arrival variance (deg 2 ) 

L = path length 

Lo = a function of density and cross-section of scattering along the 
■ path 

Measurements at The Ohio State University of the ATS-6 2, 20 and 30 GHz beacons 
as the satellite underwent synchronous orbit transition were used to derive' 
empirical constants for. this model with an effective turbulence height, h, of 
6 km (typical). r e = mean earth radius = 6377 km, 9 = elevation angle, -and L = 
path length. 

L = h 2 + 2r e h + (rgSin©) 2 ]' 5 - r e sin0. 

The constants were 


Lo = 180 km 

ai 2 = 2.6 x 10" 7 f(GHz) ?/i2 l(km) ll/6 
c 2 2 = 5.67 x 10" 6 L(km) 1 ’ 56 d(m)“ l/s 
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A plot of the variance measurement, S 2 , expressed in dB,' is shown in 
Figure 2 for four representative frequencies for a 4.6 m diameter aperture. 

S 2 is plotted as a function of elevation angle and equivalent path length 
for a 6 km high region of turbulence. 

. Figure 2 represents the average S 2 as derived from the O.S.U. 

empirical constants. However, since both ai 2 , and oz 2 may be represented in 
closed form as a function of C n 2 (Tatarski-1961) , instantaneous, diurnal, 
or seasonal values for S 2 may be found from this model given an estimate of 
the appropriate C n 2 . 

2.2.1 Applicability of the Model 

The empirical constants which were found from observed data are 

applicable for the prediction of average turbulence-induced propagation 

effects in a temperate climate, during the warmer seasons of the year, aRd 

under non-precipitating clear-air conditions. It is necessary to derive 

local estimates of C 2 for the model if these conditions are not the same. 

n 

2.2.2 Distribution of Amplitude Variance 

It is known that peak-to-peak variations of 30 N-units in the 
refractive index are expected on a time scale of days and hours (Theobold- 
1978). Corresponding fluctuations in received signal amplitude variance 
expressed in dB would be expected to be about 20 dB peak-to-peak for a 
fluctuation of 30 N-units out of an average of 345. Figure 3 shows a repre- 
sentative case of average amplitude variance at 30 GHz for a 4.6 m diameter 
aperture as a function of elevation angle. Curves for plus or minus 10 dB 

variation in C 2 about the average are shown for comparison. 

n — v 

A more exact representation of the expected distribution of ampli- 
tude variance may be obtained given measured statistics, of variance vari- 
ability about the average. Figures 4 and 5 present probability distribution 
functions of variance differences for 2 and 30 GHz earth-space signals mea- 
sured over a period of 26 days. The satellite was undergoing transition in 
elevation from 0.38° to 45° and the mean variance was removed as a function 
of elevation angle. The 90% confidence limits of 14.6 and 14.7 dB, respec- 
tively, are in good agreement with the statistics of expected refractive 
index variation. 
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Figure 2. Amplitude Variance for a 4.6m 
Diameter Aperture for 1 to 100 GHz 









AMPLITUDE VAR. (dB) 


EQUIV. PATH (KM) 


200 50 30 20 15 10 



EL. ANGLE (DEG) 


Figure 3. Effect of 20 dB Peak-to-Peak 
Variation of C n 2 on Amplitude Variance 
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Figure 4. Distribution of Amplitude Variance 
from that Predicted from Average Turbulence- 
Induced Fluctuation Theory 
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Figure 5, Distribution of Amplitude Variance 
from that Predicted from Average Turbulence- 
Induced Fluctuation Theory 




2.2.3 Power Spectral Density 

The formulation of the structure of the power spectral .density of 
turbulence-induced amplitude fluctuations has been derived from classical 
turbulence theory (Tatarski-1961 ) . The theoretical spectrum of amplitude 
fluctuations in a medium characterized by a real refractive index is found 
to roll off as f~8/3 s or -26.6 dB/decade, in frequency f. -This behavior is 
not a function of operating frequency, as long as the wavelength is small 
or on the order of the smallest refractive inhomogeneities. Deviation from 
this slope will occur due to non-stationarity of the scintillation process. 

The spectral slope was calculated for time records of 102.4 seconds 
at 2 and 30 GHz on the ATS-6 CW beacons as the satellite moved in elevation 
angle from 0.38 to 25 degrees (Baxter and Hodge — 1978). Spectral slope 
was found to be essentially independent of equivalent path length and mea- 
sured statistics were well centered about the theoretical value of -26 dB/ 
decade. Figures 6 and 7 present the probability distribution functions of 
the 2 and 30 GHz spectral slopes, respectively. Figure 8 presents the worst- 
case confidence limits of distribution of spectral slope from an average 
-26.6 dB/decade, for 50% and 90% of total time. Such an estimate may be 
used to directly find the expected fading rates and spectral components due 
to turbulence-induced amplitude scintillation. The data represents clear 
air statistics over a period of 26 days. 

2.3 LOW ANGLE SCINTILLATIONS/FADING 

Measurement of low elevation angle scintillations or fading have 
been made while ATS-6 was descending and ascending for US observers. In 
addition, measurements of non-synchronous satellites' have been made as their 
path sweeps through the atmosphere (Crane - 1977). These latter measurements 
are not reported here. 

As ATS-6 was being moved to 35° E longitude, experimenters at the 
Virginia Polytechnic Institute and State University measured the 20 GHz car- 
rier amplitude during clear and rainy tropospheric conditions from 9° ele- 
vation angle down (Stutzman, et al - 1975). During clear weather the 
average signal decreased due to clear air attenuation, but 2-3 dB scintilla- 
tions with a period of 4 minutes were observed at 9° elevation angle. As the 
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2-10 • 


145 






FREQUENCY 
RELATIVE TO AN 
ARBITRARY FREQUENCY 
F 0 (Hz) 


Figure 8. Confidence Limits of Distribution of Spectral Slope 

from Average -26dB/decade. 
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weather became hazier, the amplitude of the scintillations remained near 
2 dB but the period. decreased to 6 seconds. As the elevation angle decreased 
to 4 degrees, signal 'fluctuations averaged 7 dB; During rainy conditions, 
fades of 15 to 20 dB were observed for 10 may'hr rain rates. These high 
values are assumed due to -the -long path length at low elevation angles. 

Ohio State University also observed the descent of ATS-6 (Hodge and 
Theobold - 1976) from 42 to 2° elevation angle. The variances of the 
scintillations as a function of elevation angle are given in Figure 9 (com- 
pare with Figure 2). Because the distribution suggested a cosecant behavior, 
a minimum mean-square-error curve was fit to the data as noted in Figure 9. 

During the return of ATS-6 to 140° W longitude (ascending for U.S. 
observers) Ohio State and COMSAT conducted elevation-angle measurements over 
land, while the University of Texas observed over water. The ascending OSU 
data (Hodge, et al , - 1977) was taken at 2.075 and 30 GHz from -0.7 to 43.9° 
elevation angles (not corrected for refraction) with no discernable precip- 
itation along the path. The resulting variance is .similar to Figure 9 
except that the 2 GHz variance is about 10 dB below the 30 GHz data. OSU 
also observed that the received power level was significantly below that 
predicted by simple atmospheric path loss calculations. The power spectra 
of the scintillations, other than the general reduction of amplitude of all 
components, was similar for all elevation angles and for the presence or lack 
of cumulus clouds. 

COMSAT Laboratories (INTELSAT memo) conducted 13.2 GHz measurements 
at COMSAT and Westinghouse's Facility at Friendship Airport near Baltimore 
as the elevation angle increased from 3.6 to 14 degrees. The results from 
the Clarksburg site are given in Figure 10. The crosses denote the median 
signal levels and the minimum and maximum signal levels recorded during each 
45-minute period are indicated by the vertical solid lines. The median sig- 
nal strength was on the order' of "3 dB lower than the predicted clear sky 
signal strengths. Typical fluctuation periods were 2 seconds with up to 
15 db amplitude fluctuations. 
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Figure 10. Median Signal Level (X) and Fluctuation (I) Versus 
.Elevation Angle (Clarksburg, MD). 




The University of Texas (Vogel, et al - 1977} installed a 30 GHz 
receiver at Port Aransas, Texas, to obtain a propagation path entirely over 
water. The beacon was monitored for at least one hour each day for 16 days 
in September 1976 while the elevation angle varied from 1.5 to 17.3°. The 
results presented here are based on 3.36 hours of data (15% of the data base). 
The cumulative distributions of the log attenuation relative to the day's 
mean attenuation for an arbitrarily chosen record per elevation angle are 
shown in Figure 11. Negative values of the ordinate represent signal levels 
above the mean. A rain event was observed while the satellite was near 8. 5° 
elevation angle. Because the rain was over water, no rate was measured. 
However, the mean attenuation increased by 16 dB and the standard deviation 
increased to 2.7 dB during the rain. 

Measurements made at low elevation angles in the Canadian artic at 
4 and 6 GHz also demonstrate the magnitude of the effect at elevation angles 
near 1 degree (Strickland, et al - 1977) including the effects of refraction. 
The effects would be greater in a temperate region because of the higher 
values of refracti vity, although the measurements were made during July when 
the water vapor in the atmosphere is greatest. The cumulative statistics 
for low angle fading are shown in Figure 12. The 4GHz data was nearly 'ident- 
ical to the 5 GHz results. Table 1 demonstrates the 6 GHz link margins for 
low angle fading at 1 degree elevation angle. 

The fade durations at 6 GHz for link margins from 0 to 21 dB below 
the- long term medium are shown in Figure 13. 

The cumulative distribution for the rate of change of signal ampli- 
tude between 0.4 second samples was found to be identical for positive and 
negative-going signals. The distribution is shown in Figure 14, where the 
signal exceeded 1 dB/sec for 1% of the time and 4 dB/sec for 0.1% of the time. 

Scintillations of - 1 dB have been recorded as a cloud passes 
through the ATS-6 20-GHz beam (Ippolito - 1975). The fluctuation rate was 
fairly constant at about 16 cycles per minute. The scintillations lasted 
for about 200 seconds. In addition, short bursts of scintillations have 
been observed for 10 to 20 second intervals on clear days. Their amplitude 
reaches 3 dB pk-pk with a fluctuation rate of 20 Hz. 
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Figure 11. Cumulative Distribution of Signal Attenuation 












Figure 13. Cumulative Distributions of 6 GHz Fade Durations at 
Eureka for Various Fade Margins Relative to the Long Term Median 

Signal Level, July 1974. 
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TABLE I 


6 GHz LINK MARGINS FOR TROPOSPHERIC FADING AT EUREKA, CANADA 
ELEVATION ANGLE * 1 DEGREE 


Time Duration 

Reliability 

90 % 

99 % 

99.9 % 

Worst two- hours 

8.0 dB 

18.0 dB 

28.0 dB (Rayleigh) 

Worst summer day 

6.8 dB 

15.5 dB 

24.5 dB 

Worst summer week (5 days) . . . 

5.4 dB 

13.0 dB 

22.0 dB 

Worst month 5 July (15 days). . . 

3.8 dB 

10.8 dB 

20.3 dB 
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RATE OF CHANGE OF SIGNAL LEVEL (dB/s) 






3.0 PHASE AND ANGLE-OF-ARRIVAL VARIATIONS 


Phase variations arise due to the variable delay as the wave passes 

‘ *1 

through a medium with variable refractivity and also due to wavefront ripple 
introduced by the "lumpy" medium. The former is termed phase fluctuations, 
while the latter effect is called phase scintillations. 

3.1 ESTIMATION OF PHASE FLUCTUATIONS ON EARTH-SPACE PATHS 

The model of Muchmore and Wheel on (1955) may be used to estimate 
the r.m.s. phase fluctuation of a wave propagating on an earth-space path 
if the equivalent path length L expression of the previous section is employed. 
Figure 15 is representative of such an estimate. It assumes that the tur- 
bulent medium is contained in a height of 5 km and presents r.m.s. phase 
fluctuation for 3 and 10 GHz. 

3.2 ESTIMATE OF PHASE' SCINTILLATIONS ON EARTH-SPACE PATHS 

Accompanying the amplitude scintillations of a plane wave propagating 
through tropospheric turbulence are phase variations. According to .the 
theory of Tatars ki (1961) the mean-square phase variation over a distance 
transverse to the propagation path is: 

b*(p) ; k* u>"» 

where A is wavelength, L is the propagation path . 1 erigth through the region 
of turbulence, and C no is the surface structure constant. The constant ‘K is 
equal to 2.91 for the exponential C* .model (Tatarski - 1961) and equal to 
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4.57 from Ohio data (Theobold and Hodge-1978). This expression may be used 
to estimate the expected mean-square phase variation between two points, 
given art estimate of C no . 

An expression for .the rms phase fluctuation for a finite circular 
aperture antenna of diameter d, 





2L l AN' 


2v X 10 
X 


has been presented" by Muchmore and Wheelon (1955). The derivation employs 
a ray theory approach and assumes an exponental spatial correlation for the 
turbulence scale, is in radians,' Z is the scale length of the turbulent 
eddy, L is the path length through the turbulence, X is wavelength, and AN 2 
is the mean-square fluctuation in the refractivity N. When using this 
expression, one should only assume values of % and AN 2 " such that 


5m < % AN 2 < 500m, 

3.3 ESTIMATES OF AN6LE-0F-ARRIVAL VARIATIONS ON EARTH-SPACE PATHS 

The same ray theory approach may be used to estimate angle-of- 
arrival fluctuation. The expression for r.m.s. angle-of-arrival fluctuation 
in. radians is 

2v4T l an 2 " ^ 

a = x 10" 6 

8 [ l 

where all parameters are as previously defined. A Gaussian correlation func- 
tion for the scale of turbulence was assumed and one should impose the limits 

2 x 10” Vm < AN 2 / £ <2x10’ /in 

Figure 16 is an example for this expression, within the stated range of AN 2 /£, 
for a earth-space propagation path through a turbulent region of height 5 km. 
Note that a e 2 is directly proportional to path length and independent of oper- 
ating frequency. Also, c e decreases with increasing eddy size, Z, while phase 
fluctuation increases with increasing eddy size. 
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GST 



10-3 2 5 10“ 2 2 5 10" 1 2 5 10° 


RMS ANGLE-OF-ARRtVAL FLUCTUATION (MILL1 RADIANS) 


.Figure 16. R.M.S. Angle-of-Arri val Fluctuations for an Earth-Space Path. 




4.0 PREDICTION OF RECEIVED SIGNAL GAIN DEGRADATION 


4.1 GENERAL 


Wavefront tilt or angle-of-arri val variation results in time-averaged 
gain degradation; and phase incoherence results in instantaneous gain loss, 
excluding the atmospheric gas loss. 

4.2 ! ESTIMATION OF GAIN DEGRADATION 


The model of Section 2.2 for received signal amplitude veriance has 
also been used to derive an expression for gain reduction, R, defined by 
(Theobold and Hodge - 1978) 


R = 10 logio 


<y> 


no angle fluctuations 


R = 10 log io 


! c + 1 


/ 'b 2 ) 2 

i U of m(2 rrB*7 


! c + A 


where the constants are the same as those defined for the variance expression, 
S 2 . This value for R may then be combined with atmospheric gas loss in order 
to obtain an estimate of average received signal level for an earth- space 
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| path. ' Figure 17 presents an example of predicted signal levels for 2, 7.3, 
and 30 GHz for antenna beamwidths of 1.8°, 0.3°, and 0.15° respectively. 

Also included are measured signal levels, relative to zenith, from the ATS-6 
2 and 30 GHz (Devasirvatham and Hodge-1977) transmissions and TACSATCOM 
7.3 GHz (McCormick and Maynard-1972) beacons as the satellites were moving 
in elevation'angle. 
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REC. SJG. (dB) 


EQU1V. PATH (KM) 


200 50 30 20 15 10 



ELEV. ANGLE (DEG) 


Figure 17. Predicted and Measured Signal Level 
as a Function of Elevation Angle 
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5.0 SCINTILLATION FADING AND GAIN DEGRADATION 
DESIGN INFORMATION 


5,1 FADE DISTRIBUTION FUNCTION ESTIMATION 

The estimates of gain reduction and signal variance parameters, 

n 

R and S , have been presented. These quantities may be incorporated into 

distribution functions which are of the form used in link design. They 

represent the long term average fade statistics due to clear air amplitude 

2 

and angle-of-arrival fluctuations. The estimates of R and S may be more 
closely matched to local and seasonal conditions if a local estimate of 

O 

C n is available. A hypothetical low elevation angle fade distribution 
is presented in Figure. 18. The abscissa is referenced to the signal level 
received in the absence of turbulence, i.e., including free space loss 
and gaseous absorption. The point at which the signal level is R dB is 
also the mean of the received signal; thus, one point on the fade distri- 
bution is established. The fade distribution for turbulence induced 

j 

fluctuations is assumed to be log-normal, with mean and median being equal. 
The fade distributions resulting from the Ohio State University ATS-6 
30 GHz beacon measurements (Devasirvatham and Hodge-1977) indicate that 
this log-normal assumption is valid for elevation angles above approximately 
2°. A similar observation was made concerning the 7.3 GHz fade distribution 
above 4° elevation angle observed by McCormick and Maynard (1972). 


x 
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PERCENT TIME ABSCISSA EXCEEDED 



j SIGNAL LEVEL (TYPICAL) 

REFERENCED TO SIGNAL LEVEL RECEIVED IN' 
ABSENCE OF TURBULENCE, i.e., INCLUDING 
FREE SPACE LOSS AND GASEOUS ABSORPTION. 


Figure 18. Hypothetical Fade Distribution Function 



A fade distribution may now be produced using this assumption of 
linearity. Referring to Figure 18, it was noted that the point at which the 
received signal level is R dB represents the mean signal level. For a normal 
distribution, the mean is plotted at the 502 time abscissa exceeded point, 
indicated by (T) in the figure. One standard deviation to the right of the 
mean on a normal distribution occurs at the 15.9% time abscissa exceeded level. 
It may be easily shown that the standard deviation of received signal level, 
expressed in dB and denoted c Vd g, may be written in terms of the signal vari- 
ances S 2 . THis point, Oy dB to the right of R is denoted by @ in the figure. 
A straight line drawn between points (T) and (zj now approximately represents 
the fade distribution, referenced to the mean signal level in the absence of 
turbulence induced fluctuations. This distribution was based on small fluctua- 
tion arguments and should be employed as a lower bound when estimating a parti- 
cular fade distribution. 

Deviation of this fade distribution from the expected linear form will 
occur at small time percentages. Additional fading due to precipitation, ab- 
normal refraction, or inversions in the atmosphere will cause greater fade 
depths for the small time percentages. However, the turbulence effects, which 
are always present, are still dominant for larger time percentages. For high 
elevation angles, ' i .e. , short path lengths, S 2 will be very small .and the line 
drawn through points (T) and (IT) will be virtually vertical. 

However, the precipitation effects at the lower percentages will still 
be present for short path length cases and will become the dominant feature of 
the fade distribution. 

5.2 DOMAINS IN WHICH GAIN DEGRADATION SHOULD BE CONSIDERED 
5.2.1 Estimation of Domains '• 

The effects of amplitude and angle-of-arri val fluctuation are, of 
course, most prominent for very long path lengths and/or very narrow beam- 
widths. One may estimate whether or not gain degradation ne'ed be considered 
in a path design if elevation angle (or equivalent path length) and antenna 
beamwidth are known. Figure 19 presents regimes of average gain degradation 
between 0.5 and 3 dB and where they must be considered as a function of 
elevation angle and antenna beamwidth. 
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ANTENNA BEAMWIDTH (DEG) 


EQUIV. PATH LENGTH (KM) 



Figure 19. Gain Degradation Regimes as a Function 
of Bearnwidth and Elevation Angle 


5 7 4 


16G 





'Realized gain, that is, expected gain less gain degradation, is 
plotted as a function of antenna beamwidth (for any frequency) or equiva- 
lent aperture diameter at 30 GHz in Figure 20. All equivalent aperture 
diameters are presented for an antenna aperture efficiency of 0.6. The 
curve representing zero path length L is simply the common gain approxima- 
tion G = 41253/B 2 , where B is in degrees. Realized gain curves for path 
lengths of 50 to 300 km are plotted using the model of Section 4. Equiva- 
lent earth-space path elevation angles assuming a 6 km high homogeneous 
atmosphere are presented in parentheses. 

Notice that gain degradation due to turbulence-induced fluctuation 
is negligible for beamwidths wider than about 0.7° for all path lengths. 
Degradation effects then gradually increase as beamwidth narrows from 0.7 C 
to 0'.05° and at any particular beamwidth are approximately directly propor- ~ 
tional, in dB, to path length. As beamwidth narrows beyond 0.05°, a sat- 
uration effect occurs and the degradation becomes constant for any one path 
length. 

All design figures of Section 5 represent estimates for clear air 

effects in a temperate climate during daytime and in the .warmer months of 

the year. If a local value of C n 2 is known, more accurate values of R and 

S 2 may be obtained. If local statistics of C 2 are known, statistics of R 

n 

and S 2 may be obtained. 

5.2,2 Spatial Diversity 

Paths operating at very low elevation angles with narrow beamwidth 
antennas may experience unacceptable fading due to scintillation and multi- 
path effects. The required reliability may be regained by the use of spaced 
site diversity. A site separation greater than 300 m transverse to the 
propagation path has been suggested (CCIR Doc. 5/25, 1974-78) as necessary to 
alleviate severe turbulence-induced effects. In effect, separation on the 
order of or larger than the scale size of the largest inhomogeneities in 
refractive index along the propagation path, and especially near the surface 
where refraction is greatest, results in decorrelation. of the instantaneous 
signal fluctuations and hence improved performance. 
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The aperture effects of large antennas may be circumvented if 
several phase-locked antennas, each with- relatively wide beamwidth, are 
employed in an array to achieve the desired system gain. Of course, over- 
all fade margins will be on the order of that for a single element, but 
angle-of-arri val effects are eliminated. In addition, such an array allevi- 
ates the need to mechanically track a geosynchronous satellite, as is 
necessary with. large aperture, narrow beamwidth antennas.- 
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6.0 EXAMPLE COMPUTATIONS OF SIGNAL FLUCTUATIONS AND 
LOW- ANGLE FADING ON EARTH- SPACE PATHS 


In this section examples of the parameters described in this report 
are worked' out for a hypothetical ground station located at Columbus, Ohio 
with a 4.6m (15 ft) diameter parabolic antenna observing a 28.56 GHz 
COMSTAR beacon at 10 degrees elevation angle. Actually, the COMSTAR 
satellites are not at that low an angle, but in order to demonstrate the 
effects of gain reduction this value has been arbitrarily selected. 

6.1 AMPLITUDE FLUCTUATIONS 

The variance of the received signal amplitude is calculated using 
the expression in Section 2.2. The full half-power beamwidth B in degrees 
is 70 A/d.= 70 c/fd = (70) (3xl0 8 m/sec)/(28.56xl0 9 sec -1 )(4.6m) = 0.16 degrees 
The path length of the turbulence L is computed from h = 6km, r g = 6377 km 
and 6= 10 degrees using the equation 


L = 


h 2 + 


2r h 
e 


(r e sin6) 2 


- r sine 1 


= 34 km 

The other constants are: 
L q = 180 km 

= 1.18 x 10~ 3 
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2 


a 


2 


= 8.35 x 10 


I. = 0.17 

I = 0.83 

c 


and the signal variance relative to the average power is 


S 2 = 10 log 


I 0 


9.79 x 10^+0. 14396-0. 14297 
0.83 +0.14 


= 10 log (2.03X10" 3 ) 
10 


S -27 dB 


Mote that this agrees well with the results in Figure 2. Reference to 
Tatarski (1961) would have allowed evaluation in terms of C n 2 rather 
then the formulation by Theobold and Hodge (1978) utilized here. 

Reference to Figure 5 indicates that 50% of the time the S 2 would 
be between -24 and -30 dB, while 90% of the time S 2 would be between -20 
and -34 dB. 

The power spectrum density of the fluctuations decreases at 
26.5 dB/decade (see Figure 8). That is if one considers some lower 
frequency cutoff for the amplitude fluctuations (say 0.1 Hz) the fluctuation 
power at 1 Hz is on the average 26.5 dB below the value of 0,1 Hz. Only 
10% of the time will the 1 Hz fluctuation power be only 10 dB below the 0.1 Hz 
fluctuation power. Clearly most of the fluctuation power for clear air 
fluctuations is at the low .frequencies (less than 1 Hz). 

t 

6.2 PHASE AMD AMGLE-OF-ARRIVAL VARIATIONS 

i 

Phase fluctuations are estimated from the model of Muchmore and 
Wheel on (1955) presented in Section 3.1. Data quoted in Muchmore and 
Wheel on indicate typical values for £~ 60m = 200 feet and ANP= %. Thus 
TaFF = 30 meters. 
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For a finite circular antenna (see Section 3.2) of 4.6 m diameter, 
the rms phase fluctuation is 0.85 radians = 48 degrees. 

The estimate of the angle-of-arrival fluctuations in radians (see 
Section 3.3) is calculated to be 3.2 x 10~ 5 radians = 1.8 x ICf 3 degrees. 
This is a small number compared to 0.16 degree half-power bpamwidth of the 
antenna. Also note that the limits on AN 2 /t = 8.3 x 10” ? m~* are not 
exceeded. 

6.3 PREDICTION OF THE AVERAGE RECEIVED SIGNAL GAIN REDUCTION 

The average received signal reductipn is calculated using the 
same parameters required for calculation of the amplitude, fluctuations. 
Using the relation in Section 4.2 

R = 10 log ( 0.83 + 0.17(0.84)) 

1.0 


= -0. 12 dB 


Thus during clear weather this COMSTAR beacon will on the average be 0.12 dB 
below the value calculated considering clear air attenuation only. This 
same value could be estimated from Figure 20. 

The long-term average distribution in Figure 18 is now constructed 

(15.9% point) is found to be 2Q(log e) 

10 

is the standard deviation of the receiver volt- 
detector.. 


from R and S 2 . The point 2 

(l0“l s2 l /20 )= 0.39 dB which 
age taken from a square law 
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• 1.0 DEFINITIONS 


Depolarization refers to that effect wherein an earth-space link 
wave's polarization is altered in the troposphere. For linear polarized waves 
a vertically (horizontally) polarized wave will, after passing through a media, 
h^ve a horizontally (vertically) polarized component in addition to the initial wave. 
For circularly polarized waves a RHCP (LHCp) wave will develop into a elliptical wave. 
For frequency reuse systems based on polarization isolation this coupling reduces, 
isolation and increases "cross-talk." 

The measurement of depolarization by propagation researchers usually 
has been done utilizing orthogonally-polarized feeds on a single antenna while 
observing singly-polarized satellite signals.* This parameter is called the 
cross-polarization discrimination (CPD) or cross-polarization ratio (CPR) 
defined as (Ref. 1) 


CPD 


_ power -output from the co-pol ari zed port 


power output 
1 


from the cross-polarized port 


= (CPR) 


For perfect transmitting and receiving antennas and a perfect medium this iso- 
lation could become infinite, but with practical components some leakage is 
always present. Definitions and example calculations of depolarization terms 
have been well documented in a tutorial report by Stutzman (Ref. 2). 

Unfortunately, the system designer desires- the pross-polarization 
isolation (CPI) term defined as 

CPI (dB) = co-polarized signal power (dB) - cross-polarized signal power(dB) 
on the same channel 

Fortunately, for most levels of attenuation observed, CPI * CPD = (CPR)" 1 (Ref. 3) 

* Most experimental depolarization data has been obtained from the 11.7 GHz 
right-hand circularly polarized communications Technology Satellite (CT£) 
beacon and the 19.04 and 28.56 GHz linear polarized ATT COMSTAR satellite 
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2.0 SOURCES OF UEPOLARIZaTION 


Depolarization on earth-space paths has been observed due in the 
presence of 

e rain 

0 ice 

© snow 

e multipath 

© refractive effects 

These hydrometeor and scattering effects generate depolarization because of 
their asymmetry. For example, as raindrop sizes increase their shape departs 
from spherical and becomes oblate spheroid with an increasingly pronounced flat 
bottom. For large drop sizes a concave depression develops (Ref. 4). Polarized 
microwave energy scattered from these particles can easily be converted into an 
orthogonal polarization. 

2.1 RAIN DEPOLARIZATION 

2:1.1 Depolarization Versus Attenuation 

Correlation of depolarization with rain rate has not been too success- 
ful because of the many parameters required for these calculations. However, 

l 

experimentally and analytically (Ref. 5) it has been observed that depolarization 
can be related to total attenuation by the formula 

CPD = a - b logic ( AL ) 

where CPD is the cross-polarization discrimination in dB and AL £ is the total 
attenuation in dB due to rain (not including the clear ai r- attenuation) . 

At 11.7 GHz, the University of Texas at Austin (Ref. 6) has measured 
the cumulative attenuation and depolarization statistics for the CTS 
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operational periods from 18 October 1976 to 31 January 1978. Relating the sta- 
tistics plots for equal probabilities of exceedance the results 

CPD = 41 - 20.6 log 10 (AL ) 

were obtained with excellent accuracy. For example, for a 10 dB fade, the CPD is 
20 dB. This may be related to the rain rate and rain rate statistics using the 
effective path length and specific attenuation described in Section II at 11.7 GHz. 

In Ottawa, measurements during the July and August 1976 (Ref. 7) light 
rains with attenuation less than 3 dB tended to yield scattered results. Con- 
vective storms yielded a more consistent relation between CPD and AL g , namely 

CPD = 33.7 - 19.9 log 10 (AL ) 

For both Texas and Ottawa results, attenuations below 3 to 5 dB were 
believed partially related to ice clouds which resulted in significant depolari- 
zation without attenuation. This deploarization , sometimes referred to as 
"anomalous" depolarization, has frequently been observed synchronous satellites 
over the U.S. and Canada. 

At Blacksburg, Virginia, preliminary measurements of the CTS beacon have 
yielded (obtained by curve fitting Figure 9 of Ref. 1) 

CPD a 40 - 20 logio (AL ) 
for four major storms during June and July 1977. 

The CCIR has recommended the approximate relation 

CPD =30-20 log io (AL ) 
in CCIR Draft Report 215-3 (Rev. 76). 

Depolarization models are currently under development at several organizations. 
These .include: 

Virginia Polytechnic Institute and State University ( Refs. 17 and 18) 

e Bell Telephone Laboratories (Ref. 19) 

e University of Texas at Austin (Ref. 6) 

© Communications Research Centre (Ref. 5). 
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2 . 1.2 Frequency Scaling Depolarization Measurements 

The Virginia Polytechnic Institute and State Univeristy (Ref. 8) has 
also made measurements of the depolarization at 19 GHz vertical and horizontal 
and 28.56 GHz using the Comstar beacons. Their preliminary results for August 
1977 are given in Table 1. 


TABLE 1 

CROSS-POLARIZATION DISCRIMINATION VERSUS ATTENUATION 
(Least-Mean-Square Fits) 


Blacksburge, VA . . Elevation Angle = 0 


Frequency/Polarization 

CPD = a - b log l0 (AL^) 

11 GHz, RHCP (-CTS , 0 = 33°) 

19 GHz, vertical (COMSTAR, 0 = 44°) 
19 GHz, horizontal* (COMSTAR) 

28 GHz, vertical (COMSTAR), 

CPD = 44.7 - 22.6 log 10 (AL g ) 
CPD = 47 - 24.5 logic (AL g ) 
CPD = 37.1 - 20.0 logic (AL g ) 
- CPD = 39.4 - 15.4 logio (AL g ) 

* Note: The 19 GHz horizontal data was taken -in the month of 
' September 1977. All other data is in August 1977. 


The analysis of Nowland, et al (Ref. 5) may be utilized to show 
the expected frequency dependence of the coefficients a and b in 
CPD = a - b log l0 (AL e ). Using Equations 11 and 12 of Ref. 5 and many of 
the constants in the paper, the solid curve was derived in Figure 1. The 
dashed curve was derived using the effective path length L e = 12.82R" 0,3 
and the specific attenuation in Section II. Several experimental data - 
points are shown, but these do not correlate well with the theory (possibly 
because of the polarization dependence of "a"). The important results of 
these figures are that "a" increases with increasing frequency, while "b" 
appears to be relatively constant. In the relation CPD = a - b logio(AL g ) 
this would imply that the. CPD increases with increasing frequency, but 
because of th.e rapid increase in AL g with frequency, CPD will decrease 
for increasing frequency for moderate rain rates. . 

The revised CCIR values for "a" and ll b" are (CC-IR Doc. 5/276-E, 

21 Sept 1977) 
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VALUES OF "b" IN CPD = VALUES OF "a" IN CPD 

a-blog 10 (AL) a - b !og 10 (AL) 


7? 


O EXPERIMENTAL DATA 


60 


ELEVATION ANGLE = 45° 

HORIZONTAL POLARIZATION 

CANTING ANGLE = 10° 

CANTING ANGLE DISTRIBUTION 
tn = 0.84 


50 


40 



© Q DATA OF NOWLAND, et al 

X-.--.-x USING L 0 = 12.82R -0 ' 3 
n- n CCiR (PROPOSED,r = 10°) 


OVERTICAL, BLACKSBURG, VA 


ORHCP, BLACKSBURG, VA 
ORHCP, AUSTIN, TEXAS 


OVERTICAL, BLACKSBURG, VA 


O HORIZONTAL, BLACKSBURG, VA 

_J i I 


ID 20' • 30 

FREQUENCY (GHz) 


40 


(SAME AS ABOVE) 



FIGURE 1. Frequency Dependence of the .Coefficients • 
in the Cross-Polarization Discrimination Relation 
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a = 30 log 10 (f) - 40. logio (cosG > - 20 log (sin 2 t ) 
b = 20 

where f = frequency in GHz, 6 = elevation angle and r - polarization tilt 
angle with respect to horizontal (t = 45° for circular polarization). This 
approximation is believed valid for: 1 dB AL g _< 15 dB, 10 dB <_ CPD _< 40 dB, 

8 GHz f < 40 GHz, 10 deg. < t £ 45 deg. and 10 deg. < 6 < 60 deg. This 
estimate is provisional and should be used with caution. 

2.1.3* Elevation Angle Dependence of Depolarization 

The depolarization measurements presented earlier were obtained at 
elevation angles from 24.6 degrees in Ottawa (Ref. 7) to 49 at Austin, Texas 
(Ref. 6). The general dependence of CPD versus AL g on elevation angle 9 can 
be obtained from the theoretical results of Ref. 5. Note that both the coeffi- 
cient a and the total attenuation AL depend on elevation angle. 

c 
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The elevation angle dependent results of Ref. 5 for 11.0 and 
11.7 GHz (experimental data) are shown in Figure 2. Clearly the "a" 
coefficient is elevation angle dependent, however the experimental data 
does not confirm this fact. The "b" coefficient appears to be nearly 
independent of elevation angle and does appear to agree with the limited 
data base. 

2.1.4 Phase Variations During Depolarization Events 

The Blacksburg group (Ref. 1) has also made measurements of the phase 
difference between the co-pol ari zed and cross-polarized signal components. 

The phase has been observed to both decrease and increase by about 150 degrees 
for 3 dB fades and then does not change significantly for higher attenuations. 

The differential phase has also been observed to increase and then decrease 
in the same storm. The mechanism for this plateau at 150 degrees and why the 
sign changes remains unexplained. 

More recent theoretical results and comparison to data from the Communi- 
cation Technology Satellite and the Comstar beacons are presented in Ref. 20. 

An inexpensive self-adaptive cancellation system is also proposed and its predicted 
performance is presented in Ref. 20. 


2.2 ICE DEPOLARIZATION 

Ice particles well above the height of the melting layer may have 
significant cross-polarization effects even for small values of attenuation 
(typically below 3 to 5 dB at 11.7 GHz). This effect is believed to contribute 
to the poor correlation between the excess attenuation and the cross-polariza- 
tion discrimination at these low values of attenuation. 
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VALUES OF "b" IN CPD = VALUES OF "a" IN CPD 

a — b log 10 (AL) a-blog 1 Q<AL) 




FIGURE 2. Elevation Angle Dependence of the Coefficients 
in the Cross-Polarization Discrimination Relation 
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In Austin (Ref. 6} ice depolarization was associated with either 
thunderstorms during the summer months or with clouds in the presence of 
polar air masses during the winter. An example of the percentage of time 
that CPD was. less than or equal to the abscissa given that the excess attenua- 
tion was less than or equal to 1 dB is shown for the 18 month period from 
12 June 1976 to 31 January 1978 in Figure 3. This curve shows that 
45 per cent of the time that the CPD was less than or equal to 35 dB, there 
was less than 1 dB of attenuation; 24 per cent of the time that the CPD 30 dB 
the Al g 5 1 dB and 12 per cent of the time that CPD < 25 dB the AL fe < 1 dB. 

In contrast, using the rain depolarization relation; for 1 dB yields a 
CPD - 40 dB. Therefore systems requiring 30 dB or more CPD should expect a 
significant number of depolarization events due to ice. 

In England, simultaneous measurements of the 20-and 30-GHz ATS-6 
satellite signals at elevation angles of 22-23 degrees and 9.4 radar signatures 
have demonstrated that the bright band (0°C isotherm) does not contribute 
significantly to the depolarization because of its limited extent (Ref. 8). 

Also, it has been observed (Ref. 9) that at 30 GHz ice crystals yield 
a constant value (typically 90 degrees) of the relative phase angle between the 
crosspolar and copolar signals as a function of CPD as shown in Figure 4. The 
corresponding polar plot for a heavy rain event is shown in Figure 5. In this 
case the CPD was reduced by signal attenuation and the signal to noise ratio of 
the relative phase measurement decreased as the CPD decreased. This effect 
appears to increase the scatter of the phase angle with decrease in CPD. 

English investigators have also noted that rapid changes in relative 
phase and CPD are observed in thunderstorms and are associated with realignment 
of the ice crystal orientation by the electrostatic fields. In electrically- 
active thunderstorms, these electrostatic fields discharge rapidly resulting 
in rapid relative phase shifts of 180° and rapid decreases in CPD of '27 dB 
in 20 seconds (Ref. 10) have been observed at the occurrence of a lightning 
flash. An example of this is shown'in Figure 6 where the spikes- in the 
relative phase occur for increasing CPD and result in large phase changes. 
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PERCENT OF TIME 




ATS-6, 30 GHz 
ENGLAND 


80° 



270° 

CPD NOISE LEVEL 
IS 45 dB 


FIGURE 4. POLAR PLOT OF THE CROSS POLARIZATION DISCRIMINATION 
ARISING FROM AN ICE CLOUD 
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CPD NOISE LEVEL 
IS 45 dB 


FIGURE 5 


POLAR PLOT OF THE CROSS POLARIZATION DISCRIMINATION 
ARISING FROM A HEAVY RAIN EVENT 











2.3 SNOW DEPOLARIZATION 

Snow depolarization occurs during both the winter and summer months. 
During the summer months snow exists above the 0°C isotherm. During winter as 
the isotherm lowers the thickness of the snow increases and the depolarization 
due to rain decreases. 

In Canada (Ref. 11} tests using circularly polarized diversity radars at 
frequencies near 2.9 GHz (10.4 cm wavelength) and 16.7 GHz (1.8 cm wavelength) 
af an -elevation angle of 3.2 degrees have diagnosed storms during both summer 
and winter. During June snow occurred during a storm from 2.6 to 8.2 km 
altitude and yielded a differential phase shift of 0.36 deg/km at 2.9 GHz. 

Winter data taken at 16.7 GHz gave more variable results 0.16 to 1.17 deg/km 
for moderate to heavy snowstorms ranging in altitude from 70 m to 2.6 km. 

The mean value of differential phase shift was 0.69 deg/km at 16.7 GHz. 

2.4 MULTIPATH DEPOLARIZATION 

A sharp antenna cross-polarization diagram can allow an oblique indirect 
ray to produce a significant cross-polar component. This condition is usually 
associated with low -elevation angle events where the sidelobes are receiving 
energy from the oblique ray. Measurements have been made on terrestrial links 
at 11 GHz (Ref. 12) and 22 GHz (Ref. 13). 

The magnitude of the indirect signal reflected from the earth can be 
roughly estimated from the data in Refs. 14 and 15 (data taken near 3 GHz). 
Generally the ground-scattered wave amplitude is at least 20 dB below the 
direct wave. 

2.5 REFRACTIVE EFFECTS 

Variations in the radio refractivity (dielectric constant of 
tropospheric layers can cause rotation of the polarization plane of the 
rays -refracting through the layers. This condition will occur for layers 
which are, not perpendicular to the vertical plane containing the transmitter 
and receiver as described in Ref. 16 (in French). 
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3.0 DEPOLARIZATION STATISTICS - AN EXAMPLE 


In most cases depolarization arises due to two sources: . rain and 
ice. The statistics of depolarization events can be estimated using .the 
methodology shown in Figure 7. For rain depolarization the starting point is. 
the cumulative rain attenuation statistics computed or measured via techniques 
described in Section II. For example, the cumulative attenuation statistics 
for Rosman, NC while observing the 20GHz ATS-6 satellite beacon are shown in 
Figure 8. This data was obtained by the distribution extension technique 
described in Section II. 

Since the ATS-6 was nearly vertically polarized (within 20 degrees 
of verti cal ) the formul a 

CPD = 47 - 24.5 log 10 (AL e ) 

is selected from Table 1. At this point, two techniques exist for correcting 
this formula for frequency and elevation angle. If a large frequency or 
elevation angle scaling were required, the best approach would be to use the 
equations of Nowland, et. al . (Ref. 5). However, since no frequency scaling 
is required and since the elevation angle to ATS-6 from Rosman is 47 degrees 
compared to the Blacksburg, VA elevation angle to COMSTAR of 44 degrees, no 
scaling is needed. If the effects of polarization need to be considered, the 
equations of Ref. 5 or CCIR Doc, 5/276 -E must be used. 

Within the accuracy of the required statistics the relation 
CPD (dB) =47-25 log 10 (AL e ) 

will be utilized. Substituting the attenuation values for a given percentage 
of the year (from Figure 8) into this relation yields the CPD statistics for 
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FIGURE 7. TECHNIQUE FOR PREDICTION 
OF DEPOLARIZATION STATISTICS 
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FIGURE 8. ATTENUATION AND DEPOLARIZATION STATISTICS 
FOR ROSMAN, N.C. 
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rain shown in Figure 8. Based on these results the CPD will exceed 20 dB 
during 99.9% of the year. 

However, for low attenuation levels (high CPD's), ice will introduce 
depolarization. A crude estimate of this effect can be made based on the 
data in Figure 3. Namely, when AL g £ 1 dB about 30% more CPD events occur due 
to ice and rain, than just rain alone. This translates the data point in 
Figure 8 corresponding to AL g = 1 dB from about 0.18 percent of the year to 
( 1 . 3) (0. 18%) - 0.23% of the year (indicated by an 0 in Figure 8). Thus the 
CPD curve is raised up for low values of attenuation. For attenuations above 
3 dB the rain plus ice depolarization plots are essentially all rain effects. 

Note that one could have added in the effect of ice depolarization by 
translating the corresponding rain and ice CPD's. This is probably less 
accurate because high values of CPD are equipment dependent. Therefore an 
average value of 30% ice contribution was selected rather than converting 
each value of CPD using data shown in Figure 3. This correction is based on 
very preliminary results from a single geographic location and should be 
added only to show a trend, since the absolute value is probably in error. 
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1.0 INTRODUCTION 


If the dispersion of the propagating medium is sufficiently high, the' 
phase and/or amplitude of wide bandwidth. trnasmissions may be selectively 
altered. The tropospheric effects ahve been predicted and measured to be small, 
hpwever the ionospheric dispersion may be significant. Therefore, in this 
section of the report the ionospheric effects on the path are considered. 

Measurements are underway with the COMSTAR satellite at 28.56 +. 0*264 GHz 
by Bell Laboratories. However results of these measurements have not been 
published to date. The first available data was obtained from the ATS-6 
satellite and is presented here. 
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2.0 TROPOSPHERIC EFFECTS 


2.1 AMPLITUDE VARIATIONS 

2.1.1 Theoretical Results 

Theoretical estimates of the degradation of pulse shapes through 
rain have indicated that only minor effects are observed. The calculations 
(Ref.l) indicated that pulse distortion does not become significant until 
total rain attenuations of the order of 100 dB are encountered. Since current 
link margins do not allow such high attenuations, the link will fail due to 
signal attenuation before the pulse shape affects transmission. 

However, if required, calculations of the amplitude dispersion are 
possible for both clear air and rain attenuation based on the data and for- 
mulas in Section II of this report. 

For rain the frequency dependence of the specific attenuation A 

" i 

(db/km) is 

t 1 " R > w 

where, for example, for the frequency range from 8.5 to 25 GHz 
|| = 1. 02x10" 4 f 1-1,2 

- 0.11 f .- 1-0779 


8A_ S(aR p ) _ 
9f ~ 9f 
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For example, at 20 GHz and R = 25 mm/hr, 
a = 5. 93xl0“ 2 
b = 1.12 

|| = 7. 18x10" 3 

|| = -4. 36x10" 3 

A = 2.18 dB/km 
so 

|| = 0.23 dB/(km-GHz) 

For a typical effective path length L = 12.82R" 0,3 at 45 degree elevation 
angle 

3(AL e > 5A 

-ST-* L elr = !- 14 dB/GHz 

2.1.2 Experimental Results 

. The ATS- 6 beacons at 20 and 30 GHz were both capable of being 
modulated with +180, .+360 V +540 and +720 MHz sidetone signals. Typical 
selective fading events across the 1.44 GHz bands are shown in Figures 1 and 
2, respectively (Ref.2). These are four-second averages taken on day 270 of 
1974 just before the onset of a fade event (232000Z), at the- beginning of the 
fade event (232352Z), and before receiver lock was lost during the fade 
event (232428Z). Except for fade depths in excess of 20 dB, the accuracy of 
the attenuation measurements is +1 dB. These rain fade results, while 
representative of those taken at Rosman, do not appear to be sufficiently 
accurate for deep fades because the signal levels approach the noise floor of 
the receiver. For one-minute averages, no measureable selective fading was 
observed (Ref.' 2). 

The cross-correlation of 4 and 6 GHz signals due to low angle fading 
in the Canadian arctic was found' to be low (Ref. 3).' During a 2.5 hour 
period on the day when the fading was most severe the correlation coefficient 
was 0.34 since the 6 GHz signal experienced 55% more fades than the 4 GHz signal. 
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FIGURE 1. SELECTIVE FADING NEAR 20 GHz 
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FIGURE 2. SELECTIVE FADING NEAR 30 GHz 
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2.2 


PHASE VARIATIONS 


Phase measurements have not yielded significant' results above 
10 GHz. The phase coherent sidetone signals on ATS-6 showed only minor 
variations across the 1.44 GHz bandwidths. These variations were most 
evident for the shorter (one and four second) averaging periods compared to 
the one-minute period (Ref. 2). 

Phase measurements have been attempted for- the one degree elevation 
angle satellites observed from the arctic (Ref. 3). Unfortunately, no signi- 
ficant fade events occurred and no differential phase variations were 
recorded. 
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3.0 IONOSPHERIC EFFECTS 


3.1 AMPLITUDE VARIATIONS 

Ionospheric attenuation is inversely proportional to the frequency 
squared (Ref. 4) and is generally less than 0.001 dB at 15 GHz and' an elevati 
angle 0 = 90°. The variation is approximately related to cosecant 8. 

3.2 PHASE VARIATIONS 

The group delay due to the free electrons in the ionosphere is 

(Ref. 5) 

At = 40.3 ^ • 
cf 2 

where N g is the total electron content in electroris/m 2 , c = 3xl0 8 m/sec and f ’ 
is in Hertz. This delay is equivalent to a phase delay (in radians) 


so that 

A<f> = (2ir) ('40.3)N e /(cf). 

For a typical value of N g = 10 17 nr 2 the total phase delay at 11.7 GHz is 
only 7.21 radians. The frequency dependence of this is only 

3(A(j>) = __ 2tt(40.3 )N e 


- - 6. 2x10" 10 radian/Hertz 



=0.62 radian/GHz 
= 35 degrees/GHz. 

Of course for higher frequencies, the rate of change of phase with frequency 
decreases . 
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1.0 OVERVIEW AND SOURCES 


The absorbing gases and hydrometeors in the troposphere are a source 
of broadband incoherent noise. This is normally expressed as a sky temper- 
ature T s observed by the receiving system, such that the noise power is 
Pmatcc = kT c B n where k = Boltzmann's constant = 1.38 x 10- 23 joule/degree 
and B n is the noise bandwidth (equivalent to the bandwidth of a rectangular 
filter whose noise output is the same as the receivers actual bandwidth 
limiting filter). The other sources of noise are extra-terrestrial sources 
(primarily the sun) and noise entering the antenna sidelobes at low elevation 

angl es. 

The above noise sources are observed by ground station antennas. 
Satellite antennas pointed at a ground- station observe noise radiated 
by the earth and cloud tops. 

Hence precipitation in the path will both reduce the system carrier 
strength and raise the background noise level, resulting in a reduced carrier-to 
noise ratio. The increase in noise level (noise temperature) is particularly 
significant in low-noise receiver front ends (4 dB or less), where the resultant 
noise figure during a precipitation event could be increased by 2 dB or more. 
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2.0 SKY NOISE OBSERVED BY GROUND STATIONS 


2.1 TROPOSPHERIC CONTRIBUTION TO SKY NOISE 

The effective sky noise due to the troposphere is primarily dependent 
on the attenuation at the frequency of observation. For sky radiometers an 
empirical equation has been derived (Ref. 1): 

AL e (dB) = -10 log 10 PjlZJl 

\ T m 

where T m is the mean absorption temperature in degrees Kelvin of the atten- 
uating medium and AL e is the specific attenuation times the effective path 
length. The value (Ref. 1) 

T m = 1.12 (surface temperature in °K)-50°K 

has been empirically determined. For a typical T m = 275° K and as Ts increases 
toward T m , the radiometer tends to saturate. For typical radiometer (AT)rms 
values, the radiometer can only accurately measure to 15 or 20 dB attenuation 

i 

levels. 

' Inverting this relation, yields the resulting sky temperature for 
a given path attenuation AL e (ip dB) of 

} (AL e /-10)l 

Ts = T m \1 - 10 / 

This relation may now be utilized for estimation of the sky temperatures 
using the data in Section II of this report for AL e (dB). 
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2.1 .1 


Clear Air Sky Noise 


'The sky noise contributed by water vapor and oxygen may be -computed 
from the clear air attenuation. A typical set of curves (Ref. 2) as a func- 
tion of elevation angle 6 are shown in Figure 1. Here the ground station 
is assumed to be at sea level, one atmosphere (1013.6 millibars), surface 
temperature is 20°C and water vapor density is 10gm/m 3 (R.H. = 58 %). 

■ As an example, for a 20 GHz downlink at 10 degrees elevation angle, 

T s = 1 50°K. In a 1000 MHz noise bandwidth the noise power would be 2 x 10“ 7 
watts . 

2.1.2 Sky Noise Due to Rain 

The value of sky noise increases significantly during rain events. 
When the total attenuation approaches 10 to 20 dB the sky noise temperature 
is nearly T m . 

Actual values of T s may be obtained using the specific attenuation 
and effective path length techniques in Section II. 

Cumulative statistics of sky noise temperature versus climate region 
are now given for the Global model observing the ATS-6 20 GHz beacon from 
Rosman, N.C. The elevation angle is 47 degrees. The station elevation = 880 
meters, the 0°C isotherm ~ 3700 m above the station and so the path length 
0 = 3 km. For this path length the rain rate average parameter r = 1 . 6R“ 0 * 13 
so the average rain rate along the path is 

Rave = 1.6 R°‘ 87 

where R is the one minute average surface rain rate in mm/hr. Referring to 
Table 3.5 of Section II and the sky noise temperature relation, the values 
in Table 1 are developed. These are plotted in Figure 2 along with the rain 
rate distribution. 

The sky noise temperature (see last column of Table 1) will degrade 
the overall system noise figure of the link. For example, for a receiver with 
4 dB noise figure, the resultant noise figure for the rain rate corresponding 
to 0.01% of the year will be 5.4 dB, i.e., an increase of 1.4 dB. 


2-2 


313 



ro 

i 

co 




o 


LU 
CC 
D 
I — 
< 
cc 

UJ 

cl 


LU 


LU 

00 

O 

> 

s/ 



FREQUENCY (GHz) 


FIGURE 1. SKY NOISE TEMPERATURE DUE TO CLEAR AIR 


TABLE 1 

CUMULATIVE STATISTICS OF SKY TEMPERATURE 
DUE TO RAIN FOR ROSMAN, N.C. AT 20 GHz 
l m = 275°K 


PERCENT OF 

POINT RAIN 

AVERAGE RAIN 

TOTAL RAIN 

SKY NOISE 

YEAR 

RATE VALUES 

RATE 

ATTENUATION* 

TEMPERATUREt 

0.001 . 

102 mm/ hr 

89 mm/hr 

47 dB 

275°K ' 

0.002 

86 

77 

40 

275 • 

0.005 

64 

. 60 

30 

275 

0.01 

49 

47 

23 

274 

0.02 

35 

35 

16 

269 

0.05 

22 

24 

11 

252 

0.1 

15 

17 

7 

224 

0.2 

9.5 

11.3 

4.6 

180 

0.5 

5.2 

6.7 

2.6 

123 

1.0 

3.0 

4.2 

1.5 

82 

2.0 

1.8 

2.7 

0.93 

53 


NOTES: 

* At 20 GHz the specific attenuation A = 0.06 Rave 1 " 12 dB/km and for 
Rosman, N.C. the effective path length is 5.1 km to ATS-6. 

t For a ground temperature of 1 7 °C = 63°F the Tm' 2 * 275°K, 
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SKY NOISE TEMPERATURE (oK) 



FIGURE 2. CUMULATIVE STATISTICS FOR' SKY NOISE TEMPERATURE DUE TO RAIN 


RAIN RATE (mm/hr) 







,.2.1,3 Sky Noise Due to Clouds, Foq, Sand and Dust ' 

r 

The major contributor to the sky noise temperature is the medium 
with the highest attenuation. Generally, clouds will present the highest 
attenuations. For example, for cumulus clouds with no precipitation the 
water density will be approximately 0.5 gm/m 3 . For the Rosman example de- 
scribed earlier, (20GHz) 

ALe = K c p w t c esc 6 

where t c is the thickness of the clouds (typically 2 kilometers)'. Using 
typical numbers 

AL e = (0.4 dB m 3 /gm km)(0.5 gm/m 3 )(2km) esc (47°) 

= 0.55 dB 

The corresponding sky noise contribution is then 

Ts = T mc {l-10 (0 - 56/ ' 10) }. 

The value to select for T mc is unknown, but letting T mc equal the temperature 
of the cloud should be conservative (i.e., T mc = 0°C = 273°K). The result 
^’s only 32°K. Clearly, rain represents a much more significant contributor 
to the sky noise temperature. 

2.1.4 Total Sky Noise Temperature Arising From Several Contributors . 

The sky noise temperature from several sources do not add linearly. 
Rather, the attenuation from each contributor must be added and the total 
result substituted into the sky noise versus attenuation relation. For 
example, for the Rosman ground station observing the ATS-6, the clear air 
attenuation is 1.2 dB yielding a Ts (clear air) = 66°K. From Table 1 for 
0.2% of the year (105 minutes) the rain induced sky temperature is 180°K. 

The actual sky temperature due to clear air and rain is 203 0 K which is sig- 
nificantly less than the sum of each contributor (246 0 K) . During rain con- 
ditions the cloud contributions should also be added, but these will generally 
be even a smaller contribution than the clear air attenuation. 
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2.2 EXTRATERRESTRIAL SOURCES OF SKY NOISE 

2.2.1 Solar Noise ! 

- | 

The sun is the major source of sky noise when col linear with the 
ground station - satellite path. For geosynchronous satellites this occurs 
during the equinox periods each year for a short period each day. The sky 
noise as a function of frequency has been measured {Ref. 3). The noise level 
appears to level off at -183 dBW/Hz (active sun) and -188 dBW/Hz (quiet sun) 
above 20 GHz (see Figure 3). 

Another approach is to model the sun as a 7000°K blackbody in the 
microwave region (Ref. 4). The power per unit area and Hertz is then 

2Trhf 3 2Trf 2 kT s l W \ 

Wf c 2 (e hf/kTs -l) " c 2 W\ ~z) 

where T s = 7000°K = apparent solar temperature 
k - Boltzmann's constant 

- 1.38 x 10~ 23 Watt sec/degree 
h = 6.626 x 10“ 34 Watt sec 2 
c = 3 x 10® m/sec. 

The effective sky noise temperature observed by the ground station 
is then proportional to the integral of the solid angle subtended by the sun 
over the antenna field pattern. The sun subtends approximately a % degree 
arc when- observed from the earth. Also the sun power density is attenuated 
by the troposphere, so that during rain periods the power density is reduced. 
This is the mechanism for suntracker propagation measurements. 

2.2.2 Lunar Noise t 

The moon reflects solar radio energy back to earth. It'also sub- 

/ 

tends approximately a % degree arc from the earth with an apparent temperature 
of 205 to 21 5°K depending on frequency. However, because of the phases of 
the moon and the ellipticity of its orbit, the solid angle of the moon in a 
propagation path varies throughout the month. Reference 5 has calculated 
the necessary parameters , .but as a worst case a temperature of 21 5° K is 
typical . 
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FIGURE 3 


FREQUENCY IN MEGAHERTZ 


VALUES OF NOISE FROM QUIET AND ACTIVE SUN. 
SUN FILLS ENTIRE BEAM (REF. 3). 





2.2.. 3 Radio Stars 


The strongest radio stars are ten times weaker than the lunar emis- 
sion. The strongest stars (Ref. 6) emit typically 10 23 W/m 2<, Hz in the 10 
to 100 GHz frequency range. Three of these strong sources are Cassiopeia A, 
Taurus A and Orion A. These sources are sometimes utilized for calibration 
of the ground station G/T. During the calibrations the attenuation due to 
the troposphere is usually cancelled out by comparing the sky noise on the 
star and subtracting the adjacent- (dark) sky noise. Thus the attenuation- 
induced sky noise cancels out. 
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3; 0 NOISE OBSERVED BY SATELLITE-BORNE RECEIVERS 


Receiving (uplink) antennas aboard spacecraft observe the warm earth 
or the earth's cloud cover. This represents a noise source competing with 
the uplink signal from the ground station. 

Measurements of the brightness temperature of the earth measured 
with the electrically scanned microwave radiometer aboard Nimbus-5 at 19.35 
GHz have demonstrated brightness temperature differences of 160°K when conn 
paring water ( 1 30 0 K) and land (290°K) (Ref. 7).- Clouds appear to have a 
temperature near 270°K (0°C). Because of the high percentage of clouds and 
ocean areas on the earth, a full earth coverage antenna will observe approx- 
imately 270°K. However, when observing a ground based temperate zone ground 
station a 300°K brightness temperature is not uncommon. The power density 
is then calculated from the relation kT s B n . 
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CONCLUSION 


This report represents an ongoing effort to compile in one place a 
concise, readily useable handbook for satellite system engineers. This report 
is the initial effort and will be updated and expanded in FY 1979. 

People utilizing this report are encouraged to contact the authors 
V/ith their suggestions and criticisms of both format and content. These 
comments, where appropriate, will be incorporated in later editions of the 
handbook. 
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